FACILITY FORM 602 


NASA TECHNICAL 
MEMORANDUM 


NASA TM X- 58051 
August 1970 



SMALL OSCILLATIONS OF A VISCOUS 
ISOTHERMAL ATMOSPHERE 


A Dissertation Presented to the 
Faculty of the Department of Graduate Mathematics 
of Adelphi University 
In Partial Fulfillment of 


the Requirements for the Degree 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
MANNED SPACECRAFT CENTER 
HOUSTON, TEXAS 



SMALL OSCILLATIONS OF A VISCOUS 


ISOTHERMAL ATMOSPHERE 


Robert M. Myers 
Manned Spacecraft Center 
Houston, Texas 



PREFACE 


The study of the earth’s atmosphere has attracted 
the talents of many investigators. The author hopes that 
this research is not an unworthy contribution to this 
field of scientific investigation. 

During the conduct of this research the author has 
received considerable aid from several individuals. Since 
many numerical calculations were performed, it was neces- 
sary to develop a large number of computer programs. 

During this phase of the research, Melvin J. Ardlt and 
E. 0. Smith of the Lockheed Electronics Company were most 
helpful and the author wishes to express his appreciation. 

The author has been very fortunate since an unusual 
set of ci. cumstances* has lead to the completion of this 
research. The author wishes to express his gratitude to 
Henry P. Decell of the University of Houston and Jay M. 
Lewallen of the Manned Spacecraft Center for providing 
the opportunity to complete this research. In addition. 
Dr. Decell and Dr. Lewallen have made many helpful 


The author completed this research while on active 
duty in the U.S. Army and assigned to the Theory and 
Analysis Office of the Manned Spacecraft Center in Houston, 
Texas . 

^Formerly Chief of the Theory and Analysis Office of 
the Manned Spacecraft Center. 



- 111 - 


suggestions and comments and the author is especially 
indebted to these exceptional individuals. 

The author wishes to express his gratitude to 
Dr. Michael Yanowitch o£ Adelphi University who served 
as advisor. His patience, advice, and intuitive under- 
standing have made this undertaking worthwhile. 

The author wishes to express his deepest gratitude 
to his wife and children. Their love, willingness to 
sacrifice, and tenacious spirit have made it possible to 
complete this research. 



-IV- 


ABSTRACT 

The problem of small oscillations of an isothermal 
atmosphere is investigated. The assumptions of small, 
two-dimensional waves in a viscous, compressible, 
stratified fluid with a constant dynamic viscosity 
coefficient y leads to a linear system of two second- 
order, ordinary differential equations in the vertical 
z coordinate. 

In solving the viscous problem for small y > 0 
it is found that there is an inviscid region in which the 
solution behaves like a linear combination of inviscid 
(y = 0) solutions. Several interesting cases develop 
depending on the values of the frequency a and hori- 
zontal wave number k . The most interesting case 
concerns the viscous solution for those values of o and 
k which lead to the development of inviscid solutions 
which are wavelike in z . For this case the viscous 
solution does not satisfy the radiation condition in the 
inviscid region since viscosity reflects waves in addition 
to damping wave motion for large z . Thus, the correct 
solution of the inviscid problem consists of an incident 
and a reflected wave. As y -> 0 the ratio of the ampli- 
tudes of the incident and reflected waves approaches a 
limiting value for each fixed z in the inviscid region. 
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However, the viscous solution does not approach a limiting 

value since the reflecting layer shifts to infinity 

as y 0 and thus alters the phase of the reflected wave. 

The remaining cases are investigated and several 
numerically computed solutions are determined. On the 
basis of the computations, the validity of the lineariza- 
tion is also discussed. 
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1. INTRODUCTION 

In the study of atmospheric waves, a model which is 
frequently encountered is one of waves in an ideal com- 
pressible fluid in a half-space, with a specified vertical 
temperature profile (see e.g. [1]). For the linearized 
isothermal problem the variables are separable and a 
system of two first-order, ordinary differential equations 
in the vertical z-coordinate is obtained. 

There is little difficulty in determining a funda- 
mental set of solutions since the coefficient matrix is 
constant. However, the velocity components of both solu- 
tions are not uniformly bounded for all z and, hence, 
violate the assumptions underlying the linearization. 

For the problem of forced oscillations two condi- 
tions must be imposed to specify a unique solution. One 
condition is obtained for the vertical velocity at the 
ground by assuming that the fluid maintains contact with 
the lower boundary z = 0 and another condition must be 
obtained. 

For certain values of the frequency a and hori- 
zontal wave number k , one inviscid solution has finite 
kinetic energy in an infinite column of fluid of finite 
cross section, and the other solution has infinite kinetic 
energy. For this case it is reasonable to select the 
solution with finite kinetic energy. Thus, a unique 
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solution of the inviscid problem can be obtained by 
imposing this additional requirement. The solution for 
this case is of the form 


u(x,z,t) 


"a' 

_w(x,z,t)_ 


B 


where k > 0, a > 0 and Z is real, u(x,z,t) is the hori- 
zontal velocity component, w(x,z,t) is the vertical 
velocity component, and A and B are constants. The 
solution thus propagates in a horizontal x-direction. 

If other values of a and k are considered it is 
possible to find two solutions with oblique lines of con- 
stant phase; 


u^(x,z,t) 


i 

> 

1+ 

1 

w+(x,z,t) 


j 

+1 

I 


^Zz+i (kx±Bz-at) 


where £, k, 3 and a are all positive. Both solutions 
have infinite kinetic energy. Since one of the solutions 
has upward energy flux and the other has downward energy 
propagation, it is possible to determine a unique solu- 
tion to the inviscid problem by neglecting the solution 
with downward energy propagation. This assumption seems 
reasonable since all the energy is being supplied by the 
lower boundary. However, it must be assumed that reflec- 
tion of the wave with upward energy flux is negligible. 
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This assumption (no reflection) will be called the radia- 
tion condition. 

Recently Yanowitch [2] considered a similar 
problem: two-dimensional waves in an isothermal, incom- 

pressible fluid occupying the upper half-space with a den- 
sity distribution which decreased exponentially with 
increasing altitude. By examining the viscous problem for 
a constant dynamic viscosity coefficient p , Yanowitch 
was able to introduce a reasonable requirement on the dis- 
sipation of energy due to viscosity, which he called the 
"dissipation condition." This condition requires the 
solution of the viscous problem to dissipate only a finite 
amount of energy in a column per period of oscillation of 
the lower boundary. The dissipation condition, the no- 
slip condition (u = 0 at z = 0) , and the boundary condi- 
tion on w at z = 0 are sufficient to prescribe a 
unique solution for every small y > 0 . Yanowitch then 
analyzed the viscous problem as p tended to zero. He 
was able to show that the inviscid solution with finite 
kinetic energy is approached uniformly, on an interval 

0 < A < z < B < po as p 0 ; if the parameters a and 

k yield one inviscid solution with finite kinetic energy 
and another with infinite kinetic energy. For other 
values of a and k which yield two inviscid solutions 
with oblique lines of phase propagation and both having 
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infinite kinetic energy, he was able to show that the 
reflection of waves is not always negligible. More pre- 
cisely, he was able to show that the reflection coeffi- 
cient, K (ratio of the amplitudes of the incident and 
reflected waves) , satisfied 

-it3 

rRi = ® 

2 TT 

where 3 = , H is the density scale height and £ 

is the vertical wavelength.^ The reflection in 
Yanowitch's problem was due to a relatively thin layer 
which receded to infinity as y 0 . Thus, k did not 
approach a limiting value, i.e., arg k varied as 
y -> 0 . For an incompressible isothermal model of the 
atmosphere it appears that the inclusion of viscosity not 
only damps the wave motion for large z but is also 
capable of causing reflection. Since the wave motion is 
damped for large z , Yanowitch was also able to determine 
reasonable estimates for the amplitude of the oscillation 
of the lower boundary such that the viscous solution 
remained small for all z . Thus, viscosity also provided 
a justification of the linearization. 

^Similar results were obtained by Yanowitch [3] for 
vertical oscillations of an isothermal, compressible atmos- 
phere. In addition, Lindzen [4] obtained a similar result 
for tidal waves in an isothermal atmosphere subject to 
Newtonian cooling. 
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Clearly, compressibility should be included in a 
model of the earth's atmosphere. Thus, a natural exten- 
sion of Yanowitch's approach would be to investigate two- 
dimensional, linearized wave motion in a viscous, 
isothermal, compressible, stratified fluid occupying the 
upper half-space z > 0 . 

The incompressible and compressible models differ 
in several respects. The most notable difference is that 
the incompressible model has only a low frequency (gravity 
region) range of values for a such that the inviscid 
solutions are wavelike in z , whereas the compressible 
model has both a high frequency (acoustic region) and a 
low frequency (gravity region) range separated by an 
excluded region. In addition, the compressible model has 
inviscid solutions, for certain values of o and k , 
which are eigensolutions , or free oscillations. These 
solutions are the so-called Lamb waves and the incompres- 
sible model has no such solutions. Thus, for certain 
values of a and k a resonant situation will develop 
for the compressible model, and no resonant case is pos- 
sible for the incompressible model. 

Many of the results obtained by Yanowitch for the 
incompressible fluid model are also obtained in this 
thesis for a compressible fluid model. There is one sig- 
nificant exception: the magnitude of the reflection 
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coefficient for the case of inviscid solutions which are 
wavelilce in z depends on the horizontal wave number k . 
It was found that for the range of parameters considered, 
the reflection coefficient ic^ satisfies 

, -ttB 

^rI ^ ^ ’ 

2 H 

where 3 = — ^ , H is the density scale height, and i 
is the vertical wavelength. In general, the magnitude of 
the gravity reflection coefficient more nearly equals 
e than does the acoustic reflection coefficient. If 

an error of 15 percent is tolerated then the magnitude of 
the gravity reflection coefficient can be considered 
approximately equal to e for all k in the range of 
computations. The acoustic reflection coefficient has 
rather peculiar properties for k values which correspond 
to horizontal wavelengths of 7H to 30H. For example, if a 
horizontal wavelength of 13H (about 90 km) is considered, 
then the reflection coefficient is no longer monotonically 
decreasing as 3 increases. If the horizontal wavelength 
is greater than 30H or less than 7H, then e provides 

a reasonable approximation of the magnitude for both the 
acoustic and gravity reflection coefficients. A summary 
of the results for the case of inviscid solutions which 
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are wavelike in z is provided in figures 2 through 7 in 
Section 5.2. 

In Section 2 the differential equations for linear- 
ized wave motion in a stratified compressible isothermal 
fluid are developed. A system of two second-order, ordi- 
nary differential equations in z (the altitude) are 
obtained if the fluid is assumed viscous but thermally 
nonconducting. The inviscid equations are also developed 
in Section 2 and a fundamental set of inviscid solutions 
are found. 

The remainder of this paper is concerned with an 
analysis of the viscous problem as the dynamic viscosity 
coefficient tends to zero. The most difficult mathemat- 
ical problem encountered is the so-called asymptotic con- 
tinuation problem (Section 4) . The viscous problem is 
not completely analytically tractable and thus a numer- 
ical integration of the viscous ordinary differential 
equation was performed. Since the asymptotic continuation 
problem is inherently unstable (Section 4) , a numerical 
integration is not a trivial procedure. A modification of 

2 

For the incompressible problem, which Yanowitch con- 
sidered [2], the corresponding differential equation is 
simple enough to make possible the determination of an 
integral representation of the solution. The asymptotic 
relations are then obtained from the integral representa- 
tion. A similar procedure was attempted for the compres- 
sible problem but it was unsuccessful. 



an algorithm developed by Conte [5] is contained in Sec 
tion 4. In addition, an analysis of the error is per- 
formed in Section 4 which is not included in Conte’s 
paper. In fact, the definition of error is rather nove 
A summary of all the computations and some concluding 
remarks are in Section 5. 
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2 . FORMULATION 

The complete set of equations governing two- 
dimensional flow in a viscous but thermally nonconducting 
fluid^ are : 


,lu , 9u\ 9p 

Pl^+u— +w — > + 


9x 


dz 3x 


p Au + ^ ^ (div v) 


(la) 


Pl^+u^+w^) + || + pg 


9w 


9w 


9w 
9 z 


y Aw + ^ ^ (div v) 


(lb) 


+ div (pv) 


0 


(Ic) 


and 


3t Cpp”') » 0 . 


(Id) 


where p is the density, u is the horizontal component 
of the velocity, w is the vertical component of the 


velocity, v 


, and p is the pressure; also, y 
is the ratio of the specific heats, t is time, x is 
the horizontal space coordinate, z is the vertical 


3 

For a gas, viscosity and thermal conductivity are 
related [6, pp . 47-50] and, thus, it is inconsistent to 
assume a viscous but thermally nonconducting fluid. The 
model is very simple, but it still provides a qualitative 
description of reflection. 
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space coordinate, g is the acceleration of gravity,^ 
p is the dynamic viscosity coefficient , ^ and A is the 

„ 2 g 2 

two-dimensional Laplacian operator — - + — - . Equa- 

3x 9z 

tions (la) and (lb) are the Navier-Stokes equations, equa- 
tion (Ic) is derived by assuming conservation of mass, and 
equation (Id) is the adiabatic law which neglects second- 
order dissipative terms. 

For small deviations from equilibrium the equa- 
tions (1) can be approximated by the linearized equations: 

% It * H * T 5 I ’ ^2=^’ 

1=0 If * H * 8'’ ' ^ > t2b) 

an -- 

^ + p^(div v) + w 0 ’ *^2c) 

and + c^p^Cdiv v) - gp^w = 0 , (2d) 

where the variables without subscript are the perturbed 
values and the variables with the subscript zero are the 

^The acceleration of gravity g is assumed constant 
and, also, it is assumed to be the only body force acting 
on the fluid. 

^The dynamic viscosity coefficient p is assumed 
small but constant. 
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equilibriuin values of these variables. The equilibrium 
values can be obtained from the following equations: 



dz ^*^0 

= 0 , 

(3a) 


Po 


(3b) 

and 

2 

c 

= yRTo , 

(3c) 

where 

R is the gas constant 

and the temperature 

^0 - 

assumed 

constant (isothermal) 

• 



If a scale height H 

is defined by 



H = 

RTo 

g 

C4) 

then 

1 

Pq 

1 

H ’ 

(5a) 

or 

Po(p) = 

-z/H 

p e 
g 

C5b) 

and 

2 

C 

YgH 

(5c) 


If the Brunt-Vaisala frequency is defined by 



C6) 


then 
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i 9^u 


3w 
® 3x 


3x 


(div v) 


1 


3t 




(7a) 


and 


sfw 

3t^ 


3w 
3 z 


- c 


3 z 


_ 2 ] 
(div v) + ^ (div v) 


and 


3t 




I£ the following variables are introduced 


u(x,z,t) 


(8a) 

w(x,z,t) 

= iW(z)e"<^^-^^) , 

(8b) 

P (z) 

= po<;^)/pg > 

(8c) 

z 

= z/H , 

(8d) 

k 

II 

PC 

(8e) 

a 

II 

Q 

C8f) 

e 

= P JH 

yp^H^ 

(8g) 
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into the differential equation (7) then a system of 
second-order differential equations is obtained: 


lea 

0 p 


4 

3 



“ 

1 1 


U(z) 



W(z) 



iea\ 

3 / 


k 


leo 



r 1 


U(z) 


W (z) 


p k 


o 


Y 

-k 1 


4 ,2 
lea j k 


- 1 


Y 


- 

Y 

+ ieak^ 

Y 


U(z) 

W(z) 


(9) 


d -2 

where ' = j — and p = e . The tilde has been omitted 
dz 

from the variables p, z, k and a in equation (9). No 
confusion should arise in what follows, since only the 
dimensionless variables will be considered. 

Equation (9) might arise if the periodic solution 
is desired when the lower boundary is forced to oscillate 
with frequency o . The fluid maintains contact with the 
lower boundary and, hence, must satisfy the no-slip condi- 
tion, U[0) = 0 . Without loss in generality, the oscil- 
lation of the lower boundary can be normalized so that 
WCO) = 1 . 

The problem of forced oscillations is, of course, 
artificial and the boundary conditions 


U(0) 


0 


(10a) 
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and W(0) = 1 (10b) 

are only a part of the complete mathematical formulation. 
Only the effects of the "upper boundary" are going to be 
assessed, not the effects of the oscillating lower 
boundary. For this problem it is necessary only to inves- 
tigate a plane wave incident on the upper boundary; it is 
not necessary to examine a realistic mechanism for the 
development of such a wave. Thus, the boundary conditions 
(10) only provide a normalization of the solution of the 
viscous problem and are otherwise physically meaningless. 

In addition to the boundary conditions (10) , other 
requirements must be imposed in order to insure uniqueness 
of the solution of (9) . The lower boundary is capable of 
performing only a finite amount of work per period of 
oscillation per unit area. Thus, it seems reasonable to 
require that only a finite amount of energy be dissipated 
in an infinite column of fluid of finite cross section. 
This condition will be referred to as the dissipation con- 
dition and denoted by DC . 

The local dissipation of energy depends on the 
dynamic viscosity y and the squares of the space deriva- 
tives of u and w . Thus the DC is equivalent to 
requiring 
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J j]U’(z)l^ + 1 W (z) 1 ^ + |U(z)|^ + I W(z) I ^ dz < “ 

0 ( > 

( 11 ) 


where k is assumed to be nonzero. 

Thus ^ the complete mathematical formulation of the 
viscous problem consists of a system of two second-order, 
linear, ordinary differential equations (9), boundary con- 
ditions (10) and the DC (11) . 

For p large and e small it is expected that the 
solution of the viscous problem can be approximately 
obtained by. considering e = 0 or p = 0 . If p is set 
equal to zero in (2) and relations (3) through (6) and 
(8) are used then 


U(z) 

W(z) 




U(z) 

W(z) 


( 12 ) 


is obtained. 

The solutions of (12) are easily determined since 
the coefficient matrix is constant. The first step is to 
solve for the eigenvalues of the coefficient matrix. This 
leads to the dispersion relation 


\ + a 


0 


(13a) 
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and 


2 ^ 1 1 2 
a = - k" + 

Y Y „2 


(13b) 


where A is an eigenvalue of the constant coefficient 
matrix in (12) . The roots of the dispersion relation (13) 
are 


i. 




V 


a 


and 




- a 


(14a) 

(14b) 


If a < ^ in (13), then and are real and 

X < y < X . The solutions of the inviscid equation (12) 

2 u X 

are of the form 


u. (z) 


1 




W. (z) 
j 


2 / 1 2 
a /Y - k 


- h) 



L \ 3 


X . z 
, D 


j = 1 or 2 , 


(15) 


for o^/y - ^ 0 and 


U^(z) 


(z) 


- x: 


X + — 
1 Y 


k X, + - 
1 Y 


XiZ 


(16a) 


and 


U,(z)- 


1' 

W2CZ) 


_ 0 _ 


^2^ 


(16b) 


for — - k = 0 
Y 
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The solutions with the subscript 2 in relations 
(15) and (16) have finite kinetic energy in an infinite 
column of finite cross section, since 

(2X . -1) z 


j|Uj(z)| 


(z) I ^ jp (z) 


constant x e 


(17) 


and " 2 


The solutions corresponding to the sub- 

1 


script 1 have infinite kinetic energy since ^ J 

For the case — - k"^ = 0 the solution with finite 
Y 

kinetic energy, the so-called Lamb wave, is a free oscil- 
lation. Hence, for this case a resonant situation can be 
expected to develop. 

If ct = ^ then the dispersion relation (13) yields 
only one root. For this case the solutions of (12) have 
the forni 


U^(z) 


1 

(z) 


2 / .2 
a / Y - k 


[WTi\ 



z/ 2 


(18a) 


and 


a, 1 

to 

N 

1 


1 

z/2 ^ 

o' +7 

U^(z) 

wpz) 


(2 + y) (o^/y - k^) 

tJ ~ L 

W^(z) 


(2 - Y)k(| - i) 



(18b) 


Both solutions have infinite kinetic energy in an infinite 
column of fluid of finite cross section. 
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The most interesting case occurs when the roots of 
the dispersion relation (13) are complex. This occurs 
when 


a 


2 ^ 1 1 2 ^ 
H_ - > 1 

Y Y (j2 4 


The roots of (13a) are then of the form 


J - 13 , 


2 - 13 , 


where 


3 




a - 4 . 


Thus 3 is a dimensionless vertical wave number 

2 , 5 IT? 



Figure 1. — Inviscid dispersion relation 
for Y = 1.4. 


(19) 


(20a) 

(20b) 

(20c) 
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The solutions of the inviscid differential equation (12) 
are wavelike in z if the dimensionless parameters a 
and k have values which lie in the shaded regions of 
figure 1. The upper shaded region A in figure 1 is 
referred to as the acoustic region and the lower shaded 
region G is called the gravity region. Both solutions of 

X 

(12) for a > -^ are of the form (15) , Both solutions 
have constant kinetic energy for all z and, hence, have 
infinite kinetic energy in a column of fluid. 

If the group velocities of solutions (15) for com- 
plex Aj are investigated, then it can be shown that 


U^(z) 

W^Cz) 


2 2 

has an upward energy flux for a /y ■ k >0 and 


2 2 

a' downward energy flux for a /y - k <0 (see [7]) . 


Similarly , 


U2(z) 

W2(z) 


can be shown to have a downward 


energy flux for o^/y 


> 0 and an upward energy flux 


2 2 

for a /y - k < 0 . Thus, in the acoustic region the 
vertical phase propagation and the energy flux are in the 
same direction for solutions (15) , and for the gravity 
region the vertical phase propagation and the energy flux 
are in opposite directions. 

It is worth noting that in some publications prior 


to Yanowitch’s paper [2] the upward- and downward-going 
waves are handled separately, that is, linear combinations 
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o£ the two solutions are not considered. It appears that 
the tacit assumption regarding this wave motion is that 
reflection does not occur. 

In the following sections the differential equa- 
tion (9) will be investigated subject to the boundary 
conditions (10) and the DC (11) . The solution of the vis 
cous problem will be investigated in the limit as e 0 
(or y ^ 0) and also for a value e = 10 , which is 

comparable to the value in the earth's atmosphere. The 
object is to obtain some insight into the correct formula 
tion of the inviscid problem (y = 0) . 



- 21 - 


3. REFORMULATION AND DISCUSSION OF THE VISCOUS PROBLEM 


3.1 REFORMULATION 

It will be convenient to introduce a new indepen- 
dent dimensionless variable ^ , defined by 


The differential equation (9) then takes the form 


1 0 U(0 

^2 


0 C - 


W(5) IkU - 


k 5 - j l U(C) 


W(E) 


u- - 


Y - 1 


- - C 

Y 


^ fr ^ 1 2 

— K + k 

Y 


U(0 


W(0 


where 0 = ^ 


. The boundary conditions (10) 


become 


U(C^) _ 0 

W(5^) [i 


where 


1 

ico 


(23b) 
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Thus , as e -> 0 through positive values “ along 

the ray arg ^ . Any bounded region, 

0<A<z<B<oo, is restricted to a line segment on 
the ray arg 5 = -^ and as e 0 the line segment 
shifts to infinity. Hence, an investigation of (9) as 
e->0 on 0<A<z<B<<» is equivalent to an investi- 
gation of (22) for large ^ . 

U(C)“ 

0UC?) 


If 


iCO 


W(0 

0W(O 


(24a) 


I 


10 0 0 
0 10 0 
0 0 10 
0 0 0 1 


C24b) 


0 0 0 0 

0 0 0 0 

0 0 0 0 

000 -3/4 


(24c) 



23- 


R 


4k' 

3 

0 


0 

0 

k 


0 

0 

3k‘ 


k 

'3 


(24d) 


and 


D 


- k^ 
Y 


0 

5(y - l)k 
4y 


0 

0 

•3k 

4 


k 

Y 


3a' 

4y 


k 

0 

3 

4 


(24e) 


then (I ?K) ^ (R + ^D)y(0 


(25) 


Equation (25) has regular singularities at ^ = 0 and 
^ j and an irregular singularity at ^ = °° . The vis- 
cous problem, formulated in Section 2, can be restated in 

terms of the new ^-variable. The problem is to investi- 
gate solutions of equation (25) on the ray arg g = , 

which satisfy the boundary condition (23) and the DC (11) 
as e ^ 0 . 

Near 5=0 there exist solutions of equation (25) 
which exhibit the scalar growths 5 , (^n 5)5 > 5 > and 

(iln 5)5 ^ as 5^0 (see Appendix A) . Only the 
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Ic k: 

solutions which grow like 5 and (£n 5)? satisfy the 
DC. Hence, the solution of the viscous problem is merely 
a linear combination of the two solutions which satisfy 
the DC . 


For large ^ a fundamental set of asymptotic solu- 


tions of equation (25) can be found (see Appendix B.l) 

■^1 "^2 

which exhibit the scalar growths 5 ^ 




vf 


, where and X^ are the 


5^e2^ and 

roots of the dispersion relation (13) . Due to the dif- 
ferent scalar rates of growth the asymptotic solutions are 
significant in different regions. In particular, due to 
the boundary condition (23) and the rapid growth of 


- t/s 


as 


^1 increases along arg ^ = ~j > “the solu- 
tion of equation (25) , which exhibits this approximate 
scalar growth, is significant only near z = 0 or 
? = 5i 

Thus, if £ > 0 is small, the two solutions which 
satisfy the DC must combine so as to eliminate approxi- 
mately the rapidly growing solution over most of the 
interval 0 < |^| < and arg ^ ^ ’ Except for a 

scaling constant the viscous problem can be approximately 
solved over most of the ^"interval by determining a linear 


^Thi s 


is discussed more fully in Section 3.2. 
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combination of the solutions satisfying the DC, which is 
asymptotic to a lineai' combination of the solutions exhib- 

1 4-1 C~^2 1 ^4 -20 ?Y 

iting tne scalar growths 5 ,5 , and ^ e ' . 

In solving the viscous problem four distinct 
regions develop. Near 5=0 the viscous forces dominate 
and the solution decays like 5 and (In 5)5 and, hence, 

this region is called the viscous region. For large 5 
the solution of the viscous problem is approximately a 
linear combination of inviscid solutions and, hence, this 
region is called the inviscid region. Connecting the vis- 
cous region and the inviscid region is a transition region 
where 5 varies from large to small values or equiva- 
lently the kinematic viscosity varies from small to large 
values. In addition, near 5^_ = a boundary layer 

develops. The object is to investigate the viscous 
problem in the inviscid region as e ^ 0 . 

What is required mathematically is a means of con- 
necting the viscous behavior to the inviscid behavior 
through the transition region, or a method to connect the 
expansions about the regular singularity 5=0 to the 
asymptotic expansions about the irregular singularity 
5 = ” . Due to the very complicated three-term recursion 
relation for the expansions about 5=0, it was not 
possible to attack this problem in the same manner in 
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which Yanowitch solved the incompressible problem. 
However, it is possible to solve the viscous problem 
numerically (See Section 4) . 

3.2 HEURISTIC DISCUSSION OF THE VISCOUS PROBLEM 

In this section a heuristic discussion of the vis- 
cous problem, formulated in Section 2, is carried out. 
Unless certain pathological situations develop, it is 
shown that it is reasonable to expect existence and 
uniqueness of the solution of the viscous problem for 
every sufficiently small value of e > 0 . For the par- 
ticular values of the parameters used in the calculations 
it appears that none of the pathological situations devel 
oped. In addition, the four regions (viscous, transition 
inviscid, and boundary layer) are discussed. 

In Appendix A it is shown that there exist pre- 
cisely two linearly independent solutions of the differen 
tial equation (25) which satisfy the DC and two linearly 
independent solutions which violate the DC. The two solu 
tions which satisfy the DC will be denoted by D£]^(5) and 
DC^ (g) . No singularities exist on the ray arg ^ ~ 
for 0 < |g| < “ ; hence it is possible to analytically 
continue DC^ (g) and DC^ (g) to the whole ray. 

THEOREM 1: Let the two-dimensional vectors 

dc^(g) and consist of the first and third 
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components of the four 'dimensional vectors DC^ (g) and 
respectively. If and dc and 

atpsp are linearly independent, then the viscous problem 
has one and only one solution. 

NOTE; It is already known that and DC^ (g) 

are linearly independent since these solutions of equa- 
tion (25) are linearly independent for small values of 
f Q . For linear differential equations two or more 
solutions are linearly independent on an interval which 
excludes singularities if and only if they are linearly 
independent at a single point in the interval 
[8 , Chapter 3] . 

It is not unreasonable to assume that 
dc^C?^) are linearly independent. See theorems B3 and B4 
in Appendix B.2 for a weaker hypothesis in establishing 
the existence and uniqueness of the solution of the 
viscous problem. 

PROOF: The viscous problem consists of a differen- 

tial equation (25), boundary condition (23), and the DC 
(11) . Any solution of equation (25) which satisfies the 
DC (11) must be a linear combination of DC^ (g) and DC^ (g) 
(see Appendix A) . The only question is whether there is 
one and only one linear combination of DC^ (g) and DC ^ (g) 
which satisfies the boundary condition (23) at g = g^ . 
The boundary condition (23) imposes a requirement only on 
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tlie first and third components of DC^ (C and DC^ . 
Since and dc^CS^) are linearly independent it is 

clear that there is one and only one solution of 


c^ X dc^(^^) + X 


Q.E.D. 


It is important to note that it is reasonable to 
expect the existence of a unique solution to the viscous 
problem formulated in Section 1 only if there are pre- 
cisely two linearly independent solutions of the differen- 
tial equation [25) which satisfy the DC (11) . If only 
one solution of equation (25) satisfied the DC, there 
would be little hope of also satisfying the boundary 
condition (23) . If more than two linearly independent 
solutions of (25) satisfied the DC, then there would 
probably be infinitely many solutions to the viscous 
problem . 

Theorem 1 is not very useful for obtaining addi- 
tional qualitative information for ^ large. The region 
1 5 I > 1 is important since it corresponds to the physical 


region 0 < z < ^n 


The most important analysis 


from a physical standpoint concerns the region |C| >1 • 
In order to obtain additional qualitative informa- 


tion for ? large it is necessary to make use of the 
asymptotic expansions about the irregular singularity 
C = °° , In Appendix B a fundamental set of formal 
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asymptotic expansions is developed about the irregular 
singularity C = “ • The principal results are the fol- 
lowing : 

a. The lead terms in two of the asymptotic expansions 
correspond to the two inviscid solutions, that is, 

two of the asymptotic expansions are asymptotic to the 
inviscid solutions . 

b. There are two remaining asymptotic solutions which 

i- 

4 i 2 CT T Y 

exhibit the scalar growths C e ' . 

The solutions of equation (25) which are asymptotic 
to the inviscid solutions will be denoted by INV ^ (g) and 
INV ^ (g) . The solution of equation (25) which is exponen- 
tially increasing in as | ^ [ increases will be denoted 

by BLSOL (g) . The exponentially decreasing solution will 
be denoted by TLSOL (g) . 

If y^p(?) is the solution of the viscous problem 
for some small e > 0 , then y^p (C) satisfies 

y^,p(g) = Cp X DCp(g) + c^ X DC2(g) (26) 

for some constants c^ and c^ . Since INV p (g) , I NV ^ ( g ) , 
TLSOL (g) , and BLSOL (g) grow at different asymptotic rates, 
they are linearly independent for 0 < |g| < “ and 
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arg ^ • Hence, there exist constants d^, d^, d^, 

and d, such that 
4 

Xyp^?) = X ^^(5) + d^ X IIW^CO + d3 

X TLSOL (g) + d^ X BLSOL (g) . (27) 

•7 

The boundary condition (23) restricts the norm of 
y^p(5) near z = 0 as e -> 0 , provided it is assumed 
that there exists a unique solution to the viscous problem 
for every e > 0 sufficiently small. This additional 
hypothesis is reasonable since DC^ (g) and D^^Cg) are 
linearly independent for all g such that 0 < |g| < °° 
and arg ^ ~ ~j (see theorems B3 and B4 in Appendix B, 
Section B . 2) . 

If the first and third components of y„^(g) are 

^VP 

0(1) near z = 0 for each e > 0 sufficiently small, 
then 

/ / -ReXi\ 

dp = 0(|gj y since ||l]W^(g^)|i = o(|gJ j , 

(28a) 





/ -KeA \ 

h = ) 

since 

1 llNV2(g^) 1 

1 = o(icj d 


(28b) 


y 

The only exception is the resonant Lamb-wave case. 
The details for this case are found in Appendix B, case 4, 
and in Section 5 . 
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since 


and 


since 



(28c) 


(28d) 


where and X. and X. are the roots o£ the 

1 leo 1 2 

dispersion relation (see Appendix B) . 

As |?| is decreased (z increased) from the value 
l^ii = ^ (z = 0) , the term d^x BLSOL (g) decreases very 
rapidly since 

I |d^ X BLSOL (g) I I ^ 1^1 ^ BLSOL (g^) | | 

-■^(Re/^-Re/C) • 

X e (29) 


Since Re/|Y > Re/^ if and 

arg ^ , thus, the term | | d^x BLSOL (g) | | decreases 

exponentially in / 1 ^ j as |c| is decreased or equivalently 
as z is increased. The rapid decrease of d^x BLSOL (g) is 
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even more remarkable when it is related to changes in z , 
If ARe/l" = Re/^Y ~ R.e/5 , then 


ARe/? = 


/2 


hz/2 


/2 


ea 


(30a) 


and 


-^ARev/l" 

/f 


= EXP 


2a 



ye 


(1 


- A z/ 2 N 
e ) 




2ye 


-Az 


(30b) 


Thus , fox an increase in z of about the 

term d^x BLSOL (g) is multiplied by a factor of about ^ . 
For small positive e the solution ELSOL (g) is only 
important near the boundary z = 0 , that is , the term 
d^x BLSOL (g) in equation (27) decreases so rapidly that it 
is insignificant and can be neglected outside a thin 
boundary layer. The boundary layer thickness is 0 m 
as e ^ 0 . The solution BLSOL (g) is called a boundary 
layer solution. 


Suppose the constants c. 


and 


in equa- 


tion (26) are not known, but the vectors DC^ (g) , DC^ (g) , 
INV^(g) , INV^ (g) , TLSOL (g) , and BLSOL (g) are all deter- 


mined at g 


For small positive values of e the 


term d^x BLSOL (g) is negligible at 5 ^ p" , thus, the 


constants c^ and c^ in equation (26) satisfy 
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dl X INvJr) " d2 X INvJl 


c, X DC /i) + c X DCji 
1 1 \ 1 / 2 2 \ 1 


+ X TLSOL^ij . (31 

8 ““11 

For a value of e like 10 the equation (31) 
is correct to hundreds of significant figures and, hence, 
is certainly consistent with the accuracy of the mathe- 
matical model which has been developed. 

If d^ ^ 0 in (31) , then the viscous problem can 
be solved in two steps. First solve for the constants 

^1’ ®2’ ^3’ ®4 that 


e, X DC + e, x DC^(l| + e, x INV + e, x TLSObfi- 

1 1\1/ 2 2\l 3 l\l/ 4 \l 


• iNvUr . 


and then solve for the constants d^ and d, such that 

2 4 


the vector X.„p(C) is given by 


)Lvpf^l^ = - e3lW^(?^)) + d^ X BLSOL (g^) , 


that is, satisfies the boundary condition (23) at 

^1 ~ boundary condition (10) at z = 0 . The 

solution, TLSOL (g) , is neglected in equation (33) since it 


For the earth's atmosphere, c is comparable to 
10“H [6, Appendix 1] . 
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decays exponentially fast. For e « lO"^^ equation (33) 
is correct to hundreds of significant figures. 

Thus, for large 5 the vector is simply a 

linear combination of the solutions INV ^ (^) and I NV ^ ( g ) 
which are asymptotic to the solutions of the inviscid dif- 
ferential equation (12) . This linear combination of solu- 
tions is modified in a region very close to the boundary 
z = 0 . The region where the inviscid solutions accu- 
rately approximate INV ^ (g) and INV ^ (g) and the boundary 
layer solution is negligible is the inviscid region men- 
tioned in Section 3.1. 

For values of C such that |cl <1 the solutions 
DC^CO and DC^ (Q can be obtained from the expansions 
about the regular singularity ^ = 0 (see Appendix A, 
Section A.l). Due to viscosity the solutions decay as 
does , (£n 5)C^ as 5-^0 or e and ze as 

z . The region 0 < I?] < 1 is the viscous region 

since viscous forces dominate and force the decay of the 
solution of the viscous problem. 

Connecting the viscous region to the inviscid 
region is the transition region. In the transition region 
TLSOL (g) will be significant and for this reason it will 


be called the transition layer solution. 
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Tlie viscous problem is solved for small c > 0 
when equations C32) and (33) are solved. The equa- 
tion (33) is easily solved since the solutions INV ^ (S) > 
INV ^ (g) , and ELSOL (g) are easily computed from the asymp- 
totic expansions (see Appendix B, Section B.l). However, 
equation (32) must be solved prior to equation (33) since 
the scalar e^ in equation (32) is required in equa- 
tion (33). Accurate values for. INV ^ (Q , INV ^ (g) , and 
TLSOL (g) are obtained from the asymptotic expansions (see 
Section B.l) for large ? . Thus, in order to solve equa- 
tion (32) it is necessary to continue the accurate values 
of the asymptotic expansions to small values of ^ , that 
is , to ? “ ^ ‘ This problem of continuing the accurate 
values of the asymptotic solutions will be called the 
asymptotic continuation problem. It will be investigated 
in Section 4. 
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4 . INVESTIGATION OF THE ASYMPTOTIC CONTINUATION PROBLEM 

This section was developed as a series o£ modifica- 
tions of a paper by Conte [5] . A numerical algorithm is 
developed which is specifically aimed at solving the 
asymptotic continuation problem. 

4.1 NUMERICAL ALGORITHM 

A concept which arises frequently in numerical 
analysis is stability. Often a numerical algorithm is 
considered either stable or unstable. However, certain 
problems have "inherent instability", quite independently 
of the particular numerical algorithm used [9]. Due to 
the different asymptotic groi^ths of solutions of (25) , the 
asymptotic continuation problem is inherently unstable. 
Only the solution which grows most rapidly is easily com- 
puted with small relative error. 

It is expected that the relative error for a 
standard numerical integration scheme will grow exponen- 
tially fast for the computation of any subdominant solu- 
tion in the asymptotic continuation problem. However, if 
a certain definition of error (Section 4.2) is introduced, 
then it is found that the relative error grows only alge- 
braically fast. In the remainder of this section a 
numerical algorithm is developed. An analysis of the 
algorithm follov/s in Section 4.2. 
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In Appendix B four distinct formal asymptotic solu- 


tions of (25) about the irregular singularity ^ = oo are 


developed. The formal truncated asymptotic expansions are 

denoted by y. (5) for j = 1,2,3, or 4 and 
3 I L 

/s 

L = 0,1,2,* •• . It is shov/n in Section B.2 that v CO 

X , u 

/\ 

and y _2 0 ^^^ multiples of inviscid solutions . The 


vector l 


(C) exhibits the scalar growth ^ 


1 

4g2ay Y 


as 


^ 03 


and 


exhibits the 


scalar 


growth 



In order to solve the asymptotic continuation 

problem it is necessary to determine the solutions of (25) 
1 9 

at ^ ’ which exhibit the distinct asymptotic 

growths. Thus, it is necessary to determine vectors 
(C) at ^ ~ \ such that 


0(5) ~ y. „(0 as 5 ^ ~ , arg ? = -| , (34) 

for j = 1,2,3, and 4 . Relation (34) does not uniquely 
specify all of the vectors (?) • Since any set of four 
solutions of (25) which satisfy (34) can be used in the 
asymptotic continuation problem, it is natural to consider 
families of solutions of (25) rather than uniquely 


^The value ? ^ ^ found to be a reasonable value 

at which to determine the asymptotic connection relations 
for most of the computations. 
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defined solutions. The families are chosen so that each 
member has the same asymptotic expansion on the ray^'^ 
arg ^ ■ The family of solutions of the differential 

equation (25) which are asymptotic to y . ^(C) , is denoted 
by 

THEOREM 2: If y.(C) is a member of {y.(5)} for 

3 

j = 1,2,3, and 4 , respectively, then (?) is also a 
member of (y_j(?)} if 


14(C) 

= Y4(C) , 

( 35 a) 

Il(?) = 

14(0 + cy4(C) , 

( 35 b) 

l2(^) = 

y.2^0 + d^^(?) 

( 35 c) 

and 



13(C) = 13(C) + c^y^(?) 


( 35 d) 


where c, d, c^, c^, and c^ are arbitrary complex 
scalars . 

PROOF: The vectors y^(?) and exponentially 

dominate y^(?) as ? . Hence, any multiple of y_^(?) 

can be added to y^(?) or without altering the 

asymptotic properties of y^^(?) or ' 

^°The families are easily extended to the sector 
I arg ? I < IT . However, this extension is unnecessary for 
the immediate problem of interest. 
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Similarly, L 3 C?) exponentially dominates y^(5), 
and y^(5) as ^ ^ . Hence, arbitrary multiples 

of added to L 3 C?) without 

altering the asymptotic properties of y_ 3 C 5 )> Q.E.D. 

THEOREM 3: If y^j (C) and Y ^ (?) are members of 

(y^C?)) for j = 1,2,3, and 4 , respectively, then 



X.4(?) = I4CO , 

(36a) 


I^(?) = LiC?) + o^^(?) , 

(36b) 

and 

Y_2(?) = 12*^?) + d^4(?) 

(36c) 

l3«) = 

13(5) + ^414^?^ 5 

(36d) 

for some 

scalars c, d, c^, c^, and c^ . 



PROOF: The vectors y;^C?), , 73 ( 0 , and 

are linearly independent since they have different 
asymptotic growths. Hence, any solution y(?) of the dif- 
ferential equation (25) can be represented in the form 

4 

y(0 = EbiZUO • 

j=i 3 3 

In narticular, there exist constants c such 

^ ’ JK. 


tkat 
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However , 

I HjCs) - ij,i,«)l I ' 0/ — 


IlY.CO - ?J,,C5)II = o/~TTljij,o(«) 

and 

IIZj(5) -Y.K)|| < III.C5) -£j,i,(C)ll + l|Y.(C) 

Thus 

IHj(5) - Y,(5)|| = , C38) 

\lCl V 

for j = 1,2,3 and 4 , where L is an arbitrarily large 
integer. The theorem is an immediate consequence of rela- 
tions (37) and (38). Q.E.D. 

Suppose a numerical integration of the differential 
equation (25) in the direction of decreasing |C| along the 
ray arg ? = is performed. The solution 

increases exponentially fast (in / 1 ^ | ) , y^j^(?) and 
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exhibit algebraic growth, and (O decreases exponen- 
tially fast (in /TXT) 3-S |C| decreases. Hence, it appears 
that for a numerical integration in the direction of 
decreasing XI that the transition layer solution (y^^(5)) 
is exponentially dominant. There should be little diffi- 
culty in obtaining an accurate value of the transition 
layer solution at, say, ^ i . However, it is antici- 
pated that the remaining solutions will be more difficult 
to determine accurately. 

If the initial vectors for a numerical integration 
of the differential equation (25) are determined from the 
truncated asymptotic expansions y. ,. (^} , for 
j = 1,2,3, and 4 , then an initial relative error 

^ is introduced. The initial error can be made 

to approach zero algebraically fast as the initial 
5 ^ . Due to the inherent instability of the continua- 

tion problem, an initially small error is expected to be 
magnified. The initial error may grow exponentially fast. 
Hence, it may happen that a simple numerical integration 
to ? “ will yield no significant figures in the 
determination of (f ) , even though subsequent errors due 
to numerical integration can be made arbitrarily small by 
carrying sufficient precision in the calculations. 
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In order to solve the asymptotic continuation 
problem, it is necessary only to obtain an accurate member 
in each o£ the families (y^C^)}, for j = 1,2,3, and 4 . 

An accurate determination of a member will serve 

as an error check. Thus, an attempt will be made to 
determine a member of although it is not required 

for the solution of equation (32) . 

One additional requirement must be imposed. Linear 
independence of the particular members of the families 
{yj(C)} must be required. It will be shown in the next 
section that it is possible to obtain accurate members of 
the individual families which are essentially linearly 
dependent (ill-conditioned) if only finitely many signifi- 
cant figures are retained. In order to ensure linear 
independence of the particular members of the different 
families a canonical form will be introduced. To some 
extent the canonical form is arbitrary. The two main 
properties of the canonical form are that it 

a. Ensures linear independence. 

b. Singles out unique members in each family of 
solutions . 


!l 
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In addition, as a bonus the canonical form provides an 
effective error control which overcomes the inherent 
instability of the continuation problem in many cases. 

The transition layer solution Cy_^( 0 ) is uniquely 
defined and, thus, for any value of ? it is possible to 
determine which component is greatest in modulus. This 
component will be called the maximum component of • 

If more than one component achieves the maximum, then the 
component with the lowest index is called the maximum com- 
ponent. Thus, for a specified value of ^ there will 
correspond a unique maximum component of y^(?)- 

If is a member of then consider 

IlC 5 ) = ZiC?) + ^14(0 . ( 39 ) 

Due to theorem 2 , Y^(?) is also a member of The 

vector is said to be in canonical form at 5 q if 

YiC^q) has a zero component corresponding to the maximum 
component of y.^(?Q)- A member of {y^(5)}, which is in 
canonical form at , is denoted by Simi- 
larly, a canonical form for the family intro- 

duced . 

The vector y;,2^^’^0^ ^ solution of ( 25 ), such 

that canonical form, that is, }^2(?o’’^o^ 
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has a zero component corresponding to the maximum com- 
ponent of 

THEOREM 4: For each value of such that 

I I TT 

0 ^ I?qI ^ ^0 ~ ~2 ’ canonical vector, 

ZljC5;?Q)> defines a unique member of foi" 

j = 1 or 2 . 

PROOF: Existence of at least one solution, 

ZjC5;?Q)) is assured since any member of {y^C?)) 

(j = 1 or 2) can be reduced to canonical form at any 
on the ray arg ~ 'J • Each member of {y ^ (5) } can be 
analytically continued to the entire ray, arg ^ ~ 

0 < 1 5 1 < 00 , since there are no singularities of the 
differential equation (25) on this ray. 

All that remains is to show that any two members of 
{y^j(5)} (j = 1 or 2) , which are simultaneously in canon- 
ical form at , must be identically equal. 

Suppose ) and ) are two members of 

{yj(^)>, which are in canonical form at ^ = Cq • The 
vector (? ; Cq) -Y j (? ; ?q) is equal to a multiple of 
due to theorem 3. Thus, 

= ^ 14 ^ 5 ) • ( 40 ) 

The constant c in (40) must be zero since 
y^j (^0 ; 5 q) (5q ’ ^0^ has a zero component corresponding to 
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the maximum component o£ CC) at ^ . Hence, 

The canonical form of a particular member of 
{^^(5)} at 5 (j = 1 or 2) can be considered the 

canonical form of the family of solutions at ^ = Cq , 
since any two members of {y_j(^)) yield the same canonical 
vector at 5=5^. Hence, the canonical form is a 
property of the family and not merely a property of the 
particular member which is reduced to canonical form. 

The vectors and y 2 (?;^o^ uniquely 

defined members of and respectively. For 

every Cq it is possible to determine which component of 
is greatest in modulus. This component will be 
called the maximum component of If more than 

one component achieves the maximum, then the component 
with lowest index is called the maximum component. 

Consider 
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The vector is said to be in temporary 

canonical form^^ at 5 = Cq if c is chosen so that 


I^For computational purposes the vector 
will be used and then discarded. It is in this sense that 
the vector yT(?0’^0^ temporary value. 
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has a zero component corresponding to the 
maximum component o£ Since yj^C5;?Q) and 

are unique linearly independent solutions of 
(25), the vector Y^C^qS^q) is uniquely defined and nonzero 
for each such that 0 < | < <=° and arg ~ ~j ■ 

The vector }(,j,C?q;?q) has a component which is greatest in 
modulus. The first component to achieve the maximum is 
called the maximum component. 

Now it is possible to introduce a canonical form 
for boundary layer family of solutions. Sup- 

pose is a particular member of (y^C?)} and consider 

l3(?;?o) = y3(C) + aiiC?;Co) ^ ’ f^2) 

where a, b, and c are scalars which are chosen such that 
has three components which are zero. The three 
zero components correspond to the distinct maximum com- 
ponents of y 4 (?o), scalars 

a, b, and c can be determined in essentially a back 
substitution process by first solving for c , then a , 
and finally b . 

THEOREM 5: For each value of such that 

0 < 1?q| < “ and arg ?q " > "the canonical form (42) 

defines a unique member, y^(C;?„), of {y.(5)}, 

^3 0 ^3 
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PROOF: Any member of reduced to 

canonical form (42) for any such that arg = -j 

and 0 < < “ , The only question is whether or not 

is uniquely defined. 

Let y 3 (?;?Q} and Y^C?;?q) be any two vectors which 
are both members of and in canonical form at 

5 = 5q • From theorem 3 

^ ’ ^ 43 ) 

where c^, , and are uniquely determined, or 

l3C5;?o) - ^ 

^ c^y^CC) , (44) 

where a and b in (44) are uniquely determined from 
c^ and c^ in (43). But and y 3 (?Q;?Q) have 

zero components corresponding to the maximum components of 
Z 4 C?o^’ Thus, the constants 

a, b, and c^ in (44) are all zero and this implies 

Y3(?o;?o^ - ZsC^o’^O^ = 0 • 

Thus , 


l3(UC„) 


Z3 ’ ^0^ 


Q.E.D. 
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In order to solve the asymptotic continuation 
problem it is necessary to determine members o£ the fami- 
lies for j = 1 , 2 , 5 , and 4 at ^ ^ 

theorems 2 and 3 it is permissible to use the canonical 
vectors X.jC5Q;CQ) at ” T solve equation (32). 

The asymptotic continuation problem is inherently 
unstable. However, theorems 2 and 3 imply that the solu- 
tions of the differential equation (25) , which exponen- 
tially dominate y^ (5) for decreasing 5 , do not destroy 
the calculations. For example, consider }^^(C), for 
decreasing \^\ the only solution which exponentially domi 
nates y, (O > is y.(C). If an error is initiated at some 

, which introduces a small multiple of ^^(?), then as 
1^1 is decreased the growth of y;^(C) may swamp the calcu- 
lation. However, theorems 2 and 3 imply that arbitrary 
multiples of y^(5) are acceptable, that is, do not affect 
the asymptotic continuation problem. The only difficulty 
which y.(0 creates in the computation of y, (?) is that 
the multiple of y^(?), introduced via an error, may grow 
to such proportions that y^(?) is masked by the multiple 
of y . (?) , that is, several of the significant figures of 

■^4 

the numerical approximation of y^(?) may reflect the use- 
less multiple of y (?) , which was introduced via an error 
In order to avoid this situation it is necessary to 
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control (not eliminate) the multiple of present in 

the numerical approximation of This can be effec- 

tively accomplished by reducing the numerical approxima- 
tion of to canonical form several times over the 

interval of numerical integration. 

There are essentially five steps in solving the 
continuation problem. They are the following: 

a. Initial vectors are determined from the trun- 

cated asymptotic expansions y. ^ (O at an 
appropriate initial value of ^ (see 

Appendix B, Section B.2). 

2 

b. The interval from to ^ ~ Y divided 

into N subintervals. If the variable x 
is defined by x = /? , then the N subin- 
tervals are chosen to have equal length in x . 
The integer N is chosen sufficiently large so 
that the transition layer solution does not 
grow in norm by more than a factor of 10. 

_ 

For the case k = 0.005 , the gravity wave and 
Lamb-wave numerical integrations were terminated at 
5 = -2 X lO'^i , The DC solutions were continued from 

to 5=-2x lO^i in order to determine the asymp- 
totic connection relations. 
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c. A numerical integration^^ is performed over the 

first subinterval starting at ^ the 

initial vector ^ ^ 1?2,3, and 4 . 

At the end of the subinterval the numerical 
solutions are reduced to the appropriate canon- 
ical form. Step c is repeated for each succes- 
sive subinterval. The canonical vectors 
obtained at the end of the kt/z subinterval are 
the initial vectors for the k+lst subinterval, 

d. A numerical integration^^ from ^ ^ 

^ - Y performed. The numerical solutions 
of (25) are once again reduced to canonical 
form . 

e. The boundary layer solution Cy^ 3 (?)) obtained 
at ^ ~ ^ numerically integrated from 

5 i to C = . The vector so obtained at 

is reduced to the boundary layer canonical 
form, y 3 C 5 jj?j)> by adding the proper 

The numerical step in 5 , which is used in the 
numerical integration, is uniform in /^; that is, if H 
is the uniform step in and h the step in C ? then 
h = (/C + H)2 - 5 = 2^/i H + H2 . 

^'^The numerical step in 5 which is used for this 
portion of the numerical integration is uniform in 5 . 
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multiples of the vectors y. ^(5^) for 

j = 1,2, and 4 . The canonical form of the 
/\ 

vector y. ^ (?^) is then compared with the 

^3 ^ L X 

canonical form of the vector which has been 

integrated from ^ ^ ^ ~ 

agreement, or lack of agreement, between these 

two vectors is a good error indication. In 

addition, the approximation of 

compared with y (?,-) • 

Since the numerical solutions, which approximate 
y^^(^), > must be accurate in order to 

reduce accurately to canonical form, it follows that 

step e is an error check on all of the solutions. Since 
dominates y^(5), ^4^5) for increasing 

l^j, there should be little difficulty in numerically con- 
tinuing 73(5) from 5 i to ^ . Thus, the error 

check at ? “ is a good indication of the errors 

present at ? ~ several problems where a closed- 

form solution was available, the algorithm (steps a through 
e) yielded accurate canonical vectors for small values of 
the independent variable. In addition, step e provided a 
one-significant-figure estimate of the maximum relative 
error present in the canonical vectors. This partially 
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justifies the claim that step e is an error check for all 
the numerical solutions. 

4.2 NUMERICAL ANALYSIS 

In this section the algorithm developed in Sec- 
tion 4.1 is investigated. The approximate numerical value 
of (5) , so obtained at a fixed value^^ of 5 , is shown 
to approach a vector associated with a member of {y.(?)}. 

It will be tacitly assumed that it is possible to compute 
numerical solutions of the differential equation (25) 
which are arbitrarily close to actual solutions of (25) . 

The continuation problem, thus limitec , is simply a ques- 
tion of whether the truncated asymptotic expansions, 
y. (?) , yield sufficient accuracy to determine the vec- 
tors y^ (?) at ? = i in the limit as ?^ -> °° . Since the 
initial error in y, (?) decreases algebraically fast as 
?^ -> 00 and the solutions of the diffe rential equation (25) 
grow at different exponential rates, i: is not obvious 
that the continuation problem can be s dved in this manner. 

Of course, the more difficult problem of accumu- 
lated error due to numerical integraticn cannot be 
entirely ignored. The fifth step of tl e numerical 

■ — 

For convenience the fixed value of ? is chosen to 

be 5 = 5^, 

1 
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algoritkm, developed in Section 4.1, provides a reason- 
able check of the accumulated errors. 

Consider 

where is a solution of the differential equa- 

tion (25) for all 5 • It will be assumed in all that 

7T 

follows that arg ^ . 

For an arbitrary in equation (45) the vector 

y_. (5) will not necessarily be a member of {y^.(C)} even 
though (?) is a member of {y^j(?)). It is necessary for 
the analysis to quantify or measure in some way the dis- 
parity (error) of y_. (?) with members of {y^.(?)}. Of 

course, it is possible to call | |£^| | the error in y_jg(?) 
at ? = ?^ .^^ However, this concept of error does not 
lead to a unique value since there are infinitely many 
members in the families (£j(?)}, for j = 1,2, and 3 . 

THEOREM 6: If 

E. (?;?o) = lly.-ie(5) -y.(?)ll > 

° Yj(?)e{yj(?)} ^ 

? fixed 


at 


^®The maximum norm 

^'^If j = 4 then 

^ • 


is used throughout this paper. 
|£^|1 uniquely defines the error 
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where a-nd are defined in equation (45), then 

the function exists, is continuous in ^ , and 

nonnegative for 0 < |5| < » , 

PROOF: If j = 4 , then theorem 6 is trivial 

since {y^(5)} has only one member, y^(5). Therefore, the 
limit process, implicit in the definition of E^^(C;?q), 
is trivial. Continuity of follows since y ^ ^ ( g ) 

and are continuous in g . 

In order to show that E^^(g;g^) exists for 
j = 1,2, and 3 , it is necessary to show that there 
exists a member, y^Cg) of {y.j(C)}, corresponding to each 
value of g such that ||y_. achieves a minimum. 

Let S^(g) be a set of nonnegative numbers corres- 
ponding to each value of g . A number s is included in 
S^(g) if and only if 


for some (g) (g) } . Thus, (g) is defined for all g 

such that 0 < |g| < «> since y. (g) and {y.(g)} exist for 

D ^ D 

0 < I g I < 00 . In addition, (g) is bounded below by zero 
for each value of g , Hence, (g) has a greatest lower 

bound (g.l.b.) for each value of g . The function 


E. (g;g^) exists if the g.l.b. is achieved for a member of 
je O 

(lyo). 
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Let 

s,(0 = g.l.b. S.(0 . 

The function E. exists if s.(^)eS.(5), Since s . (5) 

je 0 j j 3 

is the g.l.b. of there exists a sequence of func- 

tions y. (i) such that for a prescribed E y. (C)e{y.(5)} 
j n — j n — j 

for n = 1,2,3, ••• and ^ j ■ 

The sequence y. CO must remain bounded for fixed 

n 

5 , since for arbitrary e > 0 prescribed 
MZ-ieC?) - ^ s CO + e for n sufficiently 

large. Hence, Illj^^C?)!! ^ llYj^C?)!! + s^CO + e • Due 
to theorems 2 and 3 the vector can be expressed in 

the form 

q„(5) = qCO * , (46) 

where Xj particular members of Cy^CO) 

and {y. CO), respectively; that is, y.C5) and y, C?) do not 

K. ^ JC 

vary with the subscript n . 

NOTE: If j = 1 or 2 in equation C46) , then the 

sum over k includes k = 4 only; that is, 

c, = = c^ = 0 . If i = 3 in equation C46) , then 

the sum over k includes k = 1,2, and 4 . If j = 

equation C46) , then c, = = c^ = c^ = 0 . 

^ ^ ^ J In 2n 3n 4n 


4 in 
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The scalars in equation (46) are uniformly 

bounded in n since the vectors are uniformly 

bounded in n and }l^(5), 72^^^’ 

linearly independent. Hence, there exists a convergent 
subsequence of the scalars c^^^ , call it c^^ . Thus, 


^kn, 


is a convergent sequence of vectors and this 


sequence converges to a vector which is a member of 
{yj(C)}. Denote the limit vector by then 

is clearly nonnegative. 


Clearly, is continuous in ^ for a par- 

ticular member of (y^CC)} since (5) is continuous. 
Suppose at ^ ^ , 


-jl 


y (O is selected so that 


De 


-jl'-’l- 


and y^^ C5)e{yj (?) > • 


The vector yj^(?) exists due to the analysis which estab- 
lished the existence of ^ 1^1 ^ “ • 

For any ^ 
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since y. ■,(?,) is not necessarily a vector which minimizes 

1 ly. (?) " y- (?,) I I • Consider an arbitrary sequence v 
3 s 2 32 n 

of C values such that v ^ f as n -> ® . The cor- 

n 1 

responding sequence of nonnegative numbers 

bounded due to (47), and hence, there exists a convergent 
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lim 


Inequality (47) implies 


V 


^ "k 

’ 0) 

E . 

{ ^ > 

' De 

\ "k 


(48) 


Relation (48) follows if ^ = v in (47) , 

If only equality in (48) can hold, then theorem 6 
is established since the implication is that every conver- 
gent subsequence E. /v converges to the same limit, 

Hk Oj 

Thus, it would follow that 


or E^^(5;5^) is continuous in the variable ^ . 
Suppose 




(49) 


and let 


e 


E . 

:e 


(Ci;Co) 


IXm y, 

E 

n, T e 
k 



Now consider the 


sequence 


of vectors y , (2) 


such that 



and y, e{y,(^)}. There exists such vectors v, (f) due to 
-D 

the analysis which establish the existence of E , (£;£). 
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I£ the differential equation (25) is written in the 

form 

— = A(^)y(5) (50) 

and 1|A(?)1| = sup ] |A(Ov| I , (51) 

I lvM=l 

5 fixed 

then it is readily seen that | ]A(5) 1 1 is bounded on 
0 < 6 < 1^1 < “ ^recall arg C = -jj- Relations (50) and 
(51) imply 


11^(5)11 < lli(?o)M +/^I|A(S)M ll)l(S)M |ds| , 

^0 

(52) 

where y;(5) is a solution of (50). An immediate conse- 
quence of (52) is 

I'^^l |a(s) I I |as| 

IlLCq)!! i llzccplle ' 


Thus , 

I 



X e 
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However , as n, 
k 

^ I 1 A (S) 1 1 1 ds 1 
n 

and e since ^ and | |A(S) | | 

k 

is bounded. Thus, for sufficiently large n 

Jc 

- !^n ■ (55) 

k 

However, C53) contradicts the definition of 

Since (53) follows from relation (49), it fol- 
lows that (49) is false. Hence, only equality can hold in 
(48). Therefore, is continuous in 5 • Q.E.D. 

DEFINITION 1: The function E. (5 5 5^) is called the 

absolute error of ^^^(C) with respect to the family 

(Zj (?) > • 

It is often more useful to deal with relative 
errors rather than absolute errors. This is generally 
true when the solutions being investigated are capable of 
growing or decaying exponentially fast. Before a useful 
definition of relative error is developed, some prelimi- 
nary analysis will be performed. 

LEMMA 1; If ^j^C?) = II Zj i ^ M > subject to 
the constraints 
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a. Ij3_C?)e{yj C?) } , 

b. ^ fixed, and 

. C. e.J5;5„) ' I iKje'Si - IjiCSO I . 

then the function exists for 0 < |C| < °° . If 

f. (?) = min I 1^. (?) I 1 subject to the constraints 

j 2 j 2 

a. 

b. ? fixed, 

then the function f. (?) exists and is continuous in ? 

1 ^ 

for 0 < |?| < " . 

PROOF: Due to theorem 6, there is at least one 

vector, y . , (?) e{y . (?) } , corresponding to each ? such 
j 1 D 

that 

. IIjij.C) - y^poll . C54) 

If there are only finitely many such vectors, then 
f^^(?) is trivially constructed. If infinitely many vec- 
tors )ijj^(?) satisfy (54), then the analysis in the proof 
of theorem 6 can be repeated to establish that f^^(?) 

exists. Similarly, f. (?) exists. In addition, the 

D ^ 

analysis in the proof of theorem 6 can be repeated to 

establish the continuity of f.„(?). 

D ^ 
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Now return to the development of a definition of 
relative error. In general, relative error is defined by 
a ratio of the form 

I I ^APPROXIMATE ' ^EXACT ^ ^ 


The absolute error I |y - y I 

''^APPROXIMATE -i-EXACT ' 

definition 1. The problem is how to choose 


is defined in 


-EXACT 


from the family {y } . The following two definitions 

are possible ways of defining a unique relative error. To 
distinguish these two concepts one will be called the rela- 
tive accuracy and the other will be called the relative 
error . 


DEFINITION 2: The relative accuracy of v, (5) with 

respect to the family {y.(5)} is denoted by A. Rnd 

is defined by 


A . 


(5;q) 


E. 





( 5 ) 


(55) 


where y. (^) is defined in (45), E. (C;f„) is defined in 
theorem 6 and f^^(5) is defined in lemma 1. 

DEFINITION 3: The relative error of y. (O > with 

— 3 e 

respect to the family {y.(C)}, is denoted by R, (?;f„) and 

3 3 ^ ^ 

is defined by 
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(56) 


where is defined in (45), E^^C5;5„) is defined in 

theorem 6, and f,„(6) is defined in lemma 1. 

J ^ 

Clearly, 




(57a) 

s ince 

cpo > fjpo . 

(57b) 


The relative accuracy A,^(C;5 q) has peculiar proper- 
ties which render it unsuitable as a criterion for a 
numerical analysis. It is possible for to be 

extremely small even though y , (^) may be useless for 

3 e 

solving the continuation problem. 


EXAMPLE; Suppose (5 q) 


-4 

2 

3 

4 


and 


zpSo) = 


* . If c - 10^° and £ = 10 


10 


ll 

1 

1 

il 


in (45), then rj 10'^° and R^g(?Q;CQ) « 10^° . 

For the above example the relative accuracy is 
small since equals yi(5Q)+£Q and agrees 

xi^ith a member of {y^(C)} to 20 significant figures. 
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The relative accuracy roughly a measure o£ 

the number of significant figures of agreement between 

y^^CC^) and some mem.ber of {y^(?)} at ?q ■ On the other 

hand, the relative error is not influenced by 

the apparent agreement between some member of 

The relative accuracy of y (^) can be made 

arbitrarily small at any fixed finite value of ^ by 

simply adding a sufficiently lar^e multiple of y_^(?) to 

The absolute error is not necessarily 

small if is small. However, can only 

be made small at a fixed finite value of 5 by making 

EieC^;? 0 ^ small. The asymptotic continuation problem will 

be considered solved if the relative error R. can 

j e 0 

be made t approach zero at ^ ~ J > 3 “ 1>2,3, 

and 4 . 

Suppose y_j(5) ai"© particular members of {y_jC?)3> 
respectively, for j = 1,2,3, and 4 ; then the vector 
error, in (45), can be expressed in the form 

3=1 

where the scalars e 


are uniquely determined. 
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THEOREM 7: The relative error satisfies 


M. 


R,„(5;5„) < 


lirr) " l"3l 


(59a) 


M, 


) < 


^22(^0^ \l"il llyi( 5 o)ll 


|^^|\|ReXi-ReX2l+c 


1^2! I II 1^3! I listen) I I I ’ 


M. 


Rse^^’^o) ^ 


^3! IIZaC^o^ll'^T? 


^32 *-^ 0 ^ 


(59b) 

(59c) 


and 


(lej ||y,(C,)|| ^ Ie,| ||y,C?,)|| 
1 ^ 3 ! llZatCo)! + |ej Il5l4(e„)ll) , (SSd) 

for 1 < |c| < |5 q| , where and X^ are the roots 

of the dispersion relation (13), the scalars e^, e^, e^, 
and are defined in (58) , the scalar c in relations 

(59a) and (59b) can be set equal to zero if X^ ^ X^ and 
if “ ^2 ’ then the scalar c can be chosen 
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arbitrarily small but positive. The constants 

and are finite and independent of Cq . The functions 

f (5) are defined in lemma 1. 

D ^ 

PROOF: The proof of this theorem involves a good 

deal of heavy analysis. For this reason the main elements 
in the proof of theorem 7 will be outlined first. 

By far the most important element in establishing 
theorem 7 is the definition of the absolute error 
E. The definition of E. ) implies 

Eie(5;?o) - l^il Milfoil + |e2I 1172(011 

+ I £3 I I Il3(?) 1 1 , 

E2^C?;Co) < iej llz.i(Oll + U2I lli2(Oll 

N3I I ll3(?) M , 

EseCOfo) ^ I-3I 1113(^^11 

and 

E4^C5;5o) < UJ ll)ip?)ll + up lUpOll 

+ U3I llzjOli + U4I 1114(5)11 • 

Theorem 7 relates the error at 5 to the initial 
error at and this is where the heavy analysis arises. 

If were bounded above by some finite constant, then 
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the theorem would follow trivially since the solutions 
(5) have only a bounded growth on a bounded interval 
which excludes C = 0 . However, ICqI can be chosen arbi- 
trarily large in theorem 7. 

For particular vectors y.(Oe{y.(C)} it is possible 
to consider two cases: 


a. 1 < |5g| < M 

b. |5g| > M 


The first case is easily handled for any value of 
M < «> . The second case can be handled by the asymptotic 
expansions if M is sufficiently large. 

LEMMA 2. If C?)e{yj C5)} tor j = 1 , 2 , 3 , and 4 
and the roots, and , of the dispersion relation 

(13) are distinct, then 


llYj (?) 1 1 

-T^ 



(60a) 


for j = 1,2, and 4 and 1 < ] ? | < °° and 


I IZsC?) 1 1 




(60b) 


for 


1 < 


l?l < 


00 


(60c) 
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PROOF 


: Since (5) and £ (O are continuous for 


1 < 1?1 < “ and f (5) > 0 , it follows that 

IlYj (C) I I 

— — is uniformly bounded for j = 1,2,3, and 4 and 
1 < |5l < Mq < 0 ° (Mq held fixed). Thus, it is only nec- 
essary to examine | |y . (?) | 1 /f . „ (?) for ? large. But for 

3 3 

? large (see Appendix B.2) 


iUC ~ Cb„ . + ••05 


■X . 
3 


(j = 1 or 2) (61) 


I JO 


Jo, 3 ^ ® 


- / 

4 20 r ■ 


(62) 


and 


'-S,3 ^ ^ 


4 -20 f Y 


4 


(63) 


where b , b , c , and c are linearly indepen- 

vJ/-L \J f ^ U^«j JL^-j 

dent . 

It will be convenient for the analysis to consider 
the following vectors 

X . 

V^(?) = ? Jj?) . for j = 1 or 2 


t/i t/i 

V3(5) = IJ 5 13(5) * 5 1,(5): 
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and 


^ 4 ( 5 ) 


4 20 
e 





Due to the relations (61) , (62) , and (63) it follows 

that 


V, (O 

b . 
-Oo 

as 

^ CO ^ 

Cj = 

1 or 2 ) , 


V 3 U) 


0,3 

^ 00 

> 

and 

V 4 CO 

->• 

'0,3 

5 ^ 00 

• 

Therefore , 

1 ,( 0 , I 2 

CO, 

I 2 CO , and 

V4(5) 

are uniformly 

bounded and 

linearly 

independent for 

1 < 

1^1 < 00 . 


CASE 1: 


lly4(?)ll 

^42 


since consists of a single unique member y;^C?)* 

CASE 2: Consider | |)^^ (?) | 1 /f ^2 C?) • Lemma 1, which 

established the existence of implies the existence 

of a function c . (C) , such that 

raxn ’ 


£i,(E) 


AC5) * c_„,„(OV^C5) 


lii(5)ll l|vpe)ll 


Clearly , 
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It is necessary only to show that I I Y.i I I 

is uniformly bounded away from zero since | |Vj^(C)| | is 
uniformly bounded. The first step is to show that 
is uniformly bounded. Consider 

Hv,(oll - l|vp5)ll llvpe) - c„.„(5)vje)il 

and llvpS) + %i„(S)y^(C)ll < IIVpOll . 

I ivpS) I I 

Thus. 2 ||V^(C)1| ’ 

but |lV^C?]ll and MV^COM are uniformly bounded, con- 
tinuous, and for all finite C ^ 0 , [ |V^(C) 1 I > 0 . As 

? ^ , I IV^C?) 11 ^ 1 1 ^ 0,31 1 hence, ] |;^(5) 1 1 is 

uniformly bounded away from zero. Thus, is uni- 

formly bounded. 

There exists a constant, c^ , such that 

I lb - c c 11 = minimum > 0 (64) 

' '—0,1 0 , 3 ' ' 

since b and c are linearly independent vectors. 

— 0 f 1 0 f 3 

Now consider 1 1 (5) (5)V^ (C) 1 1 for large K , 
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Vi(5) * 


|Vi(E) - ko,i * 


X " £ 0 , 3 ’ 


*--0,1 " ^min*-^^-0, 3^ 


(65) 


Relation (65) implies 


Vi(0 


+ c . 
min 


(?)V4(5)I 


- {lIvpE) - bo.ill 

* l‘^min^E)l IIY. 4 (E) + £(, , 3 I l} • 


( 66 ) 


The first term on the right side of ( 66 ) is bounded 
aw^ay from zero as ^ due to (64) , and the remaining 

two terms tend to zero as ^->00 since V ( 5 ) ^ , , 

V,(£) -c ^ , and ic , (^) I is uniformly bounded. 

—4 —0,3 'min' ' 

Hence, relation (66) implies I I (?) 4 ^ I 

uniformly bounded away from zero, for sufficiently large 
5 and it is already known to be uniformly bounded away 
from zero for 1 < | 5 | ^ (^^0 arbitrarily large). 


CASE 3: Consider 1 | y^2 (^^ I I /^2 2 ’ Replace sub- 

script 1 with 2 in case 2. No new analysis is required. 
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CASE 4; Consider 1 (?) 1 I /f 32 ’ Lemma 1 

implies that there exist functions, c . (^) , c . (?) and 

^ Imin 2min 


c , . (?) , such that 


-1 ^ 

/? X 5 4^-20 Vy 


^32 








Thus , 




/r?T ^32 


c 4e-2a Ty 


Vfi 


l,(5) 


lt3(«nc,„,„(ovp5)*c^„._^(5)v^(5)4c^_^,je)v^(E) 


(67) 


It is necessary only to show that 

I lv3(5)-c,„,„(5)vy5)*c^„.yuvyE)*c^„,yovy5) 1 1 

is uniformly bounded away from zero on 1 < 1 ? ] < “ since 
-- - # 

5 ^e 20 TY I |y^(^) I I is uniformly bounded. In addition, 
it is necessary to consider only large ? since the lemma 
is valid for 1 < |?| < trivially. 

The first step in establishing that 
l|V 3 (?)+c^ . (?)V^ (?) + ••• I I is uniformly bounded away from 
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zero is to establish that c, . (?) . . ( 5 ), and c, . fi) 

lmx_n ’ 2min ’ 4mxn ^ 

are uniformly bounded functions. Clearly 

- •••II i I|v3(5)I1 

and 

- I|V3(E)|| < IIVjCE) t c3__^._^(E)vyE) + •••II 

Therefore , 

^ ^ ^4rain N 

< 2l|V3(0|| (68) 

and uniformly bounded on 1 < |?1 < “ . This 

implies that the functions c. . (^) (j = 1,2,4) are uni- 
^ jmin 

formly bounded. If the functions c. . (i) were not uni- 

^ jmxn 

formly bounded, then there would exist a sequence ^ 
as n “ such that 


c , . (5 ) 

]min n 


(69) 


for some j = 1,2, or 4 . 

NOTE: c, . (r) cannot become unbounded for finite 

jmxn 

5 since V^(?), Y. 4 ^^) linearly independent 
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£or all finite ? . Hence, relation (68) implies that 
c, . (E) are bounded for all finite 5 • 

jitiin 

If more than one subscript j satisfies (69) , 
consider the subscript j for which 


then 


(70) 


for infinitely many subscripts n , where k = 1,2, and 
4 , If more than one subscript j satisfies (70) , then 
choose the first subscript j which satisfies (70), Now 
consider for this subscript j a subsequence n^ such 
that 


1 c , , / 

5 \ 1 

> 



j min \ 



kmin(^ njl 

y 


for k = 1,2, and 4 . Let 


(71) 


d 


k 



then (71) and (72) imply 
there exist constants c 



such that 


(72) 


In addition. 


|Vj (”) 




k=l ,2 ,4 




= minimum > 0 , (73) 


Vj^(“°) 


lim^ rr) 


where 


(74) 
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The limit in C74) exists, due to the asymptotic 
properties o£ (5) » 3 £ j 2 , 3 , and 4 . The vectors, 

V^(oo)j V^C”), C~) , and (“) are linearly independent; 

thus, (73) follows. 


I|vj 5 


l-n 


k=l,2,4 


‘k “Jl 


1 |Vj (~) + 


E 

kf!j 


-d«n 


■ vy. 


k=l ,2 ,4 


^ E 

k=l,2,4 


■k ( ^ n , 




(75) 


However, relation (73) implies that the right side 
of (75) is positive for sufficiently large ^ ; that is, 
it is bounded away from zero for sufficiently large 5 


since | 1 V ( 5 


M’n. - 


-V 0 as n^ and 


\ / 

d < 1 . However, this implies 


.IK ' + 

3\ 








^^3 
k=l,2,4 

is uniformly bounded away from zero , or 


c . . IK 

jmin \ n 




k=l,2,4 



'“kmin 
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as n„ ^ . This contradicts ( 68 ), Hence, relation (69) 

X/ 

is false, that is c . . (^) are uniformly bounded for 

’ ^min 

j = 1,2, and 4 and 1 < | C 1 < • 

If the c, . ( 5 ) are bounded, then consider 

jmo.n 

11^3(0 * - c,„i„(OvyE) 1 1 


> llV_3(“) + c^V^(«>) + + c^V^(“)|l 

- - hwil 

♦ - ^ 2 ” II 

" l=4„i„«)l 11^(0 - I 4 WIII ■ 

The term in braces on the right side of (76) tends to zero 
as 5 . The first term on the right side of (76) is a 

positive constant. Hence, 1 I ^3 (C) (^)V^ (C) + • * * | 1 
is uniformly bounded away from zero. Q.E.D. 

LEMMA 3; If (C)e{yj (C) } for j = 1,2,3, and 4 
and the roots of the dispersion relation (13) are equal, 
that is, ^2 " 7 ’ 


llTjC?)!! < M <00 


(77a) 


for j = 1,2, and 4 and 1 < |C| < °° a-^d 



76- 




WJW' 


C77b) 


for 1 < 1 5 1 < 
PROOF : 


jLiC 5) ~ (ko,i " •••’5 


2 


i,(E) 


jljCS) 


ko,2 * I 5 b„,i' * •••! 


2 




/ rr 1 r4^20 f Y 

■=-0,3 " ^ 1 , 3 /-^" •••1« " 


y4C^) 


VI 


, r- 1^4 -2a T Y 

=0,3 " £1,3/'^* •••'5 ® 


The analysis for conclusion (77a) is the same as 
the proof for (60a) if 


ViCC) = 5 LiCO , 


V^U) 




In ^ ^ ^2 


and V 3 C?) and V^(^) are defined in the proof of lemma 2 . 
The analysis in lemma 2 will then establish conclusion 
(77a) for lemma 3. 
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In order to establish conclusion (77b) of lemma 3, 
consider 


V2(0 


j l-n 5^1 (E)) • 


The vectors ^ 4 ( 5 ) remain the same, that is, 

V^(5) and V^(C) are defined in the proof of lemma 2, No 
other modifications need be made in the analysis of 
lemma 2 to complete the proof of lemma 3. Q.E.D. 

LEMMA 4; If (?) e (C) > for j = 1,2,3, and 4 


ands and 

^2 

are the roots of 

the dispersion 

relation 

(13) , then 






ll>li(E)ll 

1 I)L2 C?) I 1 “ 

^12 

1 • 1 
TTYTTEpTT 

\TTTj 

1 ReX^-ReX^ 1 +c 

d 

CO 

1 1^2^53 1 1 

^21 

1 llzCEo) 1 1 

/I5„f 

1 ReX -ReX 1 +c 

) 

, (78b) 

riyiCrJTT - 

1 I 1 

\Trr ; 

I 


' Iy4(?) II “ I !y4(?Q)TT 


(78c) 
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and 


1 172 ^ 5 ) 1 1 

^ ^24 


II 14(C) 1 1 

TTy^TcpTT 

1 IZaCC) 1 1 

< M 3 ^ 

1 1 1 

1 IliCC) 1 1 

TTy^Tc^TTf 

1 17.3(0 1 1 

< M33 

ll!l 3 ( 5 o)ll 

TTzT^OTT 

nypqpTT 

1173(011 
TT ^ ^ N I 1 

< M . 

1 173(^0) 1 1 


TIiTUlTI 


34 


-4 '-^0' 


(78d) 


C78e) 


C78£) 


(78g) 


where the relations C78) are satisfied uniformly for 
1 < |r| < oo and the constants M. . can be determined 

13 

independently of [^^1. If ^ , then the constant 

c in relations (78a) and (78b) can be set equal to zero. 
If ^ ^2 ~ the constant c can be chosen to 

be an arbitrarily small, positive number. 

PROOF; If 1 ?q| were bounded above by some finite 
constant, then lemma 4 would be trivial. Thus, it is 
necessary to consider only |5g| large. 

For an arbitrary e > 0 prescribed, there corre- 
sponds a finite constant M(^) such that for |cj > M(e) , 


MyiC?)ll 




1 


e 


< 


< 1 + e 


(i = 1 or 2) ; (79) 
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for f 


1!liC5)1I 

1 - e < ± — < 1 + e , (X^ = X^); 


( 80 ) 


N^,i« 


1 - e < 


I I 


>’0,2 * 2 




< 1 + e 


(X, = X^); 


( 81 ) 


1 - £ < 


ly3(?) 1 1 


- 

4 20 ?• 


< 1 + e 


and 1 - e ^ 


1 ^ 0 , 3 ^ ^ 

I ly4(?) II 


(82) 


< 1 + £ 


-‘^0,35 ® 


- 

4 2a 7 Y 


(83) 




Relations (79) through (83) are simply the basic 
asymptotic properties o£ the formal truncated asymptotic 
expansions. Thus, for |^| > M(e) it is necessary to 
consider only the asymptotic estimates since 




1 + £ 


llyi(5)ll 1 . e Nii,ot^)N 
1 |y. (?) 1 1 - 1 - e I |y /o(?J II ’ 

J j f ^ 

(84) 


for 1^1 > M (c) . 



where “ and ^ ^2 ’ Relation (85) follows from 

relations (79) and (84) . Since the solutions of (25) are 
capable of only a bounded growth on a bounded interval, it 
is clear that can be chosen sufficiently large so 

that relation (85) is valid for 1 < |5| < “ . In addi- 
tion, relations (79) and (84) imply 


1 ly^C?) M 

> M 1 

nY^TflTT 

- ^2 1 


ReX^-ReX^ 


for 0 < < “ . Thus , 


lli:iC5)ll 

<; 

1 1 1 

Pi\ 


2 

TTz^imr 


mpfpTT 

'W/ 




Similarly, if “ ^2 ” T ’ then 

llz.i(«l! Ilypepll /iSolf 
llipoir - ” Ily^CCo)!!' \TeT/ 

since | ^,n ^ | <. M | ^ | ^ 


for 1 < |c| < and arbitrarily small c > 0 . Hence, 
conclusion (78a) follows. Similarly, conclusion (78b) 


follows . 
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Now consider I 1/ I M • 


and 


I |y.i (?) 1 1 

1 1 >14 (?) i ' 


< 


M,|5| 


ReXj^ + - 


,Re 2 a 



IliLiteJll 

TTYTOTf 


> 


M 4 I 5 I 




Re20 

e 



( 86 ) 


C87) 


Relations (86) and (87) are a consequence o£ (79), (83), 
and (84) . The constants and satisfy 

0 < < “ ; for the proper choice of and 

relations (86) and (87) are satisfied for 1 < |ci < “ • 
Thus , 


lliqoli 

1114 (Sill 



2a 


ReA^ + — 

" " Ke^(vr-/q> 

e 


TTy^TT^TT 


(88) 


However, for arbitrary A the function 


1^0 

,Trr 


R 0 X , + 


CT Ee^(/f-/q;) 

e ^ is uniformly bounded for 


1 < l?l < 1?, 


< 03 


since the factor 


1^0 

,Trr 


1 


R 0 A + 

I 1 4 


IS 


capable of only algebraic growth as |5| is decreased from 


y Q I and e 


Ref|(/t-/q) 


decreases exponentially fast as 



82 


1^1 is decreased from |?q|- Hence, relation (88) implies 

MyiCOll ^ llliC^o^ll 

I ly_4(?) 1 1 “ TTTTC^PTT 

for 1 < 1^1 < , or conclusion (78c) is established. 

Similarly, conclusion (78d) can be established. For the 
case ^ ^2 ~ analysis of (78d) must be modified. 

The function ly^"| — 1 e > where c is an 

arbitrarily small positive constant, must be considered. 
Since this function is uniformly bounded for 
1 < |?| < |5 q 1 < “ conclusion (78d) is established in all 
cases . 


In order to establish conclusion (78e) consider 

173 c?) I 1/ 1 ll^C?) 1 I • 


ReX 


'-3 . ,, I ■ 4„Re2afY 

ywiT - ” 5 !'’ 


4 


and 


11513(5)11 „ 2 

1 lK O il ^ "el'l « 


where 0 < Ih < < °° and 1 < E < °° . Hence , 

6 5 ' 
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I 1 1 
nyriCTT 







but 




>ReA, + 


1 


1^4 Re4^(/5-/5_) 
✓ Y 0 


is uniformlY bounded for 


1 < Id ^ Uni < 


Thus , 


llz.3(?)ll My3(?o^ll 

TTyfrorr " ^ TTz^id^irr 


for 1< 1^1 < 15 q 1 <“> , and conclusion 
lished. Similarly, conclusion (78f) can 

Consider | | CO 1 I / I 1 14 (5) 1 1 • 


ly3C^) 1 1 

’Ty^c?) I T 


M^e 


Re4a 




(78e) is estab- 
be established. 


and 


1 1 L 3 O) 1 1 Re4aVy 

TTirCOTT - 


where 0 < M_ < M„ < “ and 1 < | d < “ • Hence , 

O / 


IlljCOll Ra^(/C-/C„) 1IJL3(C„)II 

1114(5)11 - Mg ® nYirOtT 
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but e ° < 1 for 1 < U| < UqI < “ • Q.E.D. 

Theorem 7 is a consequence of lemmas 2, 3, and 4. 
Actually, estimate (59c) can be improved but it is unim- 
portant in what follows. 

Now consider the error in (45) . Due to the 

asymptotic nature of ^ ® suffi- 
ciently large L as ICqI 5 that is, if is the 

vector error associated with the initial vector y. , (5„) . 

^3 , L ^ 0 

/\ 

LEMMA 5: If y. t (?) formal truncated expan- 

sion computed by method I in Appendix B.l and y^^C^g) in 

/s 

relation (45) is equal to y. (5_) , then e in (45) 

^'2 f Li (J — U 

satisfies 



where y . (f ) e{y . (?) 1 ; that is, 
3 3 


IjC?) ~ ■ 

In addition, the e in relation (58) satisfy 

K. 


^il llyi(?o^ll 

I lYj (?o^ 1 1 



ti+ 


2 


(89a) 
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and 


I IZ^Tcp TI — 



c 




9 


£ 3 ! My. 3 (?o^N 

llyj(5o)M 

£4! 1174(^0^11 

1 |y.j (5 q^ 1 1 



(89b) 


(89c) 


(89d) 


For j = 1,2,3, and 4 . The constant c in relations 
(89a) and (89b) can be set equal to zero if the roots of 
the dispersion relation are distinct and set equal to an 
arbitrarily small positive constant if = T • 

PROOF : 

I hoi I 

llq(5o)ll “ 



as ^ 00 . Relation (90) is a consequence of the asymp- 

totic nature of the formal truncated expansion 
It will be convenient for the analysis to introduce the 
following vectors : 

X . 

V (O = ? V.(5) (j = 1 or 2, X ^ X ) , 

“J J J. ^ 

= irfr (\ = h ■ i) ■ 
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and 

that 

TTZTTs 

and 

Since 

j = 1. 

bounds 

V^(C) 

13 ( 5 ) 


V3(5) 


# 


5 4e-2a7Yy^(5) ^ 


14(0 


-i- ^ 


There exist constants c^, c^, c^, and c^ such 


= " SY.2«o’ * * SLiteo) 


1 Yj (5„) I I 




the vectors (C) are uniformly bounded for 

2,3, and 4 , the lemma will be established if 

can be found for Ic.l. 

' 3 ' 

Consider the vectors (5) defined by 




V^t?) if ^ 

1 

5^/y2(?) - I ^y^L^o) if X^ = X 


2 2 




IVjiK') + v^CO 
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and 


V4(5) 


¥3(5) - v^(C) 


then the vectors V . are uniformly bounded since 

■i( \ ^ 

~ 5 h. + •••' 


1,(0 ~ K 


1 

■2 ( 


-0,1 j 


h,2 " 2 


if A3 = A3 


1 

2 ’ 


) + • • . 

' 


•X . 


y.CO ~ ? • ^ + ••• if X^, j = 1 or 2 , 


- #( 

yOO ~ 


1 ) 

—0 , 3 — 1 , 3 j 


and y^CO ~ 5 


4-2oly )_ 


c + c — + • 
- 0,3 - 1,3 


Thus, the vectors V (^) are uniformly bounded and 


constants d, exist such that 
k 




where 

c-i = d^ , c^ = , if ^ ^2 ’ 
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d £n 5 

C2 - In ^ - 2 5 if = j\ 

S = "*4 * '^5 ; 

and c. = -d^+/^d^ 

The lemma will be established if it can be shown 

that 


I dfc 1 „ 

( l |e„l l / l liljUoJHT - ” 

for 1 < ICqI < “ . This can be established if the 
inverse of V(5), where 

v(5) = [v^(0, v^CO, 13(5), v^(o] , 

is shown to be uniformly bounded. However, V~^(C) can be 
constructed from cofactors of V(5) divided by det V(5). 
Since V(5) is uniformly bounded and continuous, it follows 
that the cofactors and the determinant of V(?) are con- 
tinuous and uniformly bounded. In addition, the determi- 
nant of VC5) is strictly nonzero since the columns of 
V (O linearly independent. For large 5 the columns 
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o£ VC 5 ) approach the vectors i -5 ’ - 5 ? 

and these vectors are linearly independent (see lemma B1 

in Appendix B.l). Hence, 


det V(5) - det 


as 5 • Thus, det V(^) is uniformly bounded away from 

zero . Q . E . D , 

THEOREM 8; The asymptotic continuation problem ca.n 
be solved, 

PROOF: Lemmas 2, 3, and 5 and theorem 7 imply that 


^le ’ ^0^ 


I ReA -Re A „ I +c 

dlEol " Me 


2 


0 




L + 1 

o' (91) 


^2e ’ ^0^ 




Re A ^ -Re A ^ +c 


u, 


V 

2 


u, 


L 1 


(92) 




1-L 

2 


(93) 


and 






L 

’2 


(94) 
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£or 1<|?|<|5 q| <“. However, L can be chosen 
arbitrarily large in relations (91) through (94) . 

Hence, Rj^(?;^Q) can be made arbitrarily small by choosing 
L and sufficiently large. Q.E.D. 

Thus, in principle, the continuation problem can be 
solved by continuing the initial vectors CC„) to 

5 i for sufficiently large L and . Therefore, it 

would appear that only a numerical integration is required 
to obtain approximate values for }^j(?) at ^ i . 
Theoretically this is true, but in practice this proce- 
dure may require more precision than is feasible. 

For the theoretical investigation of the error 
E^^(5 ;Cq), it was possible to neglect errors which propa- 
gate in an unbounded manner as Cq tends to infinity. 

For example, consider E^^(C;?q); then, due to theorems 2 
and 3, it seemed reasonable to neglect errors which ini- 
tiate multiples of y^ (5) , the transition layer solution. 
Thus, from relation (58), it is easily shown that 

• 11l 3(?)1I • However, ||y:j^(5)M and 

essentially the same rates as |^| is decreased and 

1 8 

If the roots of the dispersion relation differ by 
an integer, then the formal asymptotic solutions are not 
completely developed in Appendix B.l (see case 3 in 
Appendix B.l.) 
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11)13^5)11 decreases exponential!)^ fast. Hence, if 
1 e, I llyv(^?n)ll made sufficient!)^ small it is expected 
that can be made small. The analysis which cul- 

minated in theorem 8 justifies this heuristic argument. 

From a numerical or practical point of view the 
evvov 1 1 11 ) 1 ^( 5)11 cannot be ignored since the solution 

gtows exponentially fast as 1 5 1 is decreased. Thus, 
the term le^l 11)14(5)11 could greatly exceed 11)14(5)11 if 
I^qI is large and 1^1 << 15 q 1 • From a theoretical point 
of view it is unimportant how large |e^l | 1)14(5)1 1 
becomes, but from a numerical standpoint it is extremely 
important. Assuming that the numerical calculations are 
capable of maintaining a small relative accuracy (recall 
that relative accuracy is roughly a measure of the number 
of significant figures maintained in the calculations), 
then it is necessary to append some process which controls 
the multiple c^y;^(^) present in the approximate value of 
y^( 5 )* It is of the utmost importance to ensure that 
1^4! 1114(5)11 cannot greatly exceed 11)14(5)11 i£ only 

finitely many significant figures are maintained in the 
calculations. Otherwise, all of the significant figures 
may merely reflect the useless vector £ ^ 7 4 (5) • 

Similarly, in order to compute approximate values 
of 73(5) it is necessary to control the growth of £474(5), 
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tliat is, limit the multiple of the transition layer solu- 
tion present in the approximate solution Z. 2 eC?)- In order 
to compute approximate values of X. 3 C?) it is necessary to 
limit the multiples of y^( 5 ), present in 

This is difficult, since each of these solutions 
dominates for |^1 decreasing and intermediate 

rounding and truncation errors initiate multiples of 
y^( 5 ), canonical form which was 

introduced in Section 4.1 (see theorems 4 and 5) provides 
the necessary error control. 

If 

- y.(5)M < lO”''! ly, (O I I (96) 

J ^ J J 

where y. (C) is defined in (45) , then 
3 ® 

A.^(?;?q) < 10“^ . (97) 


Relation (97) would require at least N significant 
figures to be carried in the calculations. Even if (97) 
can be maintained, no useful bound can be placed on 
RjeC5;5o) since 


R . 

je 


(? ; 5q) 


< 


lO”** 


I Iq(E) 


(98) 
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where f , C5) is defined in lemma 1. Only for the case 
3 

j = 4 , the transition layer solution, can [97) be used 
to establish that 

< 10“^ (99) 

since {}[^C^)} consists of a single unique member }[^(^). 

For j = 1,2, or 3 in relation (98), no useful bound 
results since the families {y)j(?)} contain arbitrarily 
large multiples of }(.(?)■ Hence, | | y^, (O 1 | /f . „ (5) is in 
general unbounded, that is, not uniformly bounded for the 
family. Since it is desired to make R^^(5;5 q) small for 
the numerical calculations, it is necessary to append some 
condition which ensures that | | (?) | | /f (O is bounded 
uniformly on 1 y | ? | < °° . 

THEOREM 9: If the roots of the dispersion relation 

are distinct, then 

I |y.C5;0| I 

^ < M (100) 

f.yc) 

for j = 1,2, and 3 and 1 < |?| < °° where yj(?;C) a^-e 
the canonical vectors (see theorems 4 and 5) and the func- 
tions f . (?) are defined in lemma 1. 

J ^ 

PROOF: Let Y(?) be a fundamental solution of the 

differential equation (25) of the form 

Y(?) = 


[Yi (?) , I 2 ’ ^3 ’ ^4 ^ 
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where CP C5) > • If v ^ (?) > and 

VC?) = [v^C?), )^3(?). )UC5)] , 

then = }l 4 (?) » 

Li(?) - )Li(^) = ^i)C4(^) ’ 

- v^(?) = ’ 

and X. 3 CO - v^CC) = a)^^(5) + by^C?) + cy^^ (O • 

Thus, det Y(0 = det V(5) 

if V^(C) is defined by 

= [Yi(?;?), Yt(5;?)> Y3(?;?), i4(?)l 

and 

where y.(C;^) is the canonical form of y^. (?) for 
3 1 

3 = 1 or 3 and y;^(^;C) is the temporary canonical form 
defined in relation (41). The vector y^(^;C) differs from 
y ^2 (C) only in that multiples of y^^(?) and y^^(?) are added 
to to generate y_^(g;5). The vectors Ij^^^^C?) 

satisfy 
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and for each fixed 5 , 

for some y . f?) e{y ■ C?) ) • However, y, . (^) is not included 

in since 


for one particular member (?) C?) } . Only for 
j = 4 does 


Lemma 1 assures the existence of y . , (?) . In addition, 

d. t Y (?) E det Y . (?) = det Y(?) . 

c min 


RECALL : 


is the maximum norm. 


^Tn -i n ^ ^ ^ ^ - ^ ■ I I ^Imin I I I 1 y^ni-i n (O II 


min 


2m in 


-3min 


X y (^) 1 

^4min ' 


|det Y^(5)l 


llli(5;5)ll llrTU;«ll Ill3(5;5)ll 

>= 1114(5)11 
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since det is the product of only the maximum com- 
ponents of Y_rpil',0, and Thus, 


1 l7j_(?;?) 1 1 

^ V _ 

1 17^(5;?) 1 1 

X - 

1 173(?;?) 1 1 

iiwTOT 

r 

niiiiKiT 

r 

M>^3minLOll 


X 


I IL 4 CC) M 
TTuhinTT 


< 4! 


(101) 


The first, third, and fourth factors in rela- 
tion (101) are all greater than or equal to one. Thus, 
relation (101) will provide a bound for each of these 
terms if I 1 1 I / I I ^2min I I bounded away from zero 

uniformly for 1 < lc| < • 

Consider the function f^(?) defined by 


min 

f (?) = overall I 1 7 . (?) +c (?) y (? ) +c (Q y (? ) 1 [ . 

^ c^(?) and c^C?) 


It can be shown (analysis similar to theorem 6) 
that f^(?) exists and is continuous in 5 £ 01 ' 

1 < I ? 1 < . In addition. 


5 ^ii(0 ^ ko,i 

E y re) b 

^ - 0,2 






4 2afY- 




- 0,3 


and 
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as ^ -y CO (see Appendix B.2) and ^ i j ^ 3-nd c 
are linearly independent. By analysis similar to lemma 2 
it follows that 


)l2(Oll < M£^(?) , 


but 




Hence , 




or 


- M 




(1021 


Thus, from relations (101) and (.102) 


ipcic) 




< 4 ! X M 


for j = 1 or 3 and 




< 4! 


Clearly, the subscripts j = 1 and j = 2 can be inter- 


changed when the roots, and , of the dispersion 

relation are distinct. This amounts merely to a renum- 
bering of the inviscid solutions. Therefore, if 
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I 

’T7 


-Imin 
long as X ^ X 


is uniformly bounded for arbitrary X^ as 

I 172(^5 I I 

then TTrr r- r is uniformly 




bounded for arbitrary ^2 ^ ^ 


If A. 


1 

2 


1 

then 


\z^a;i') 1 1 


Q.E.D . 


M 

£n ? I +1 


C103) 


for 1 < ? I < 


Thus , 


for 


lliiye;s)ll 


T (|ta 5| + 1) 


1 or 3 and 


lliiye;;)|| 

Mil2mln(^El 


c C I £n ? I + 1) 


C104a) 


(104b) 


Since | S,n grows very slowly there is very little addi- 
tional difficulty for the case ^ ^2 ” T ’ 

Thus, in all cases the canonical form limits the 
ratio I ly^ (?;?) I 1/ I I 1 . In so doing, if 

has a small relative accuracy, then it also has a small 
relative error. 

The solutions of the differential equation (25) 
grow by only a bounded factor on any bounded interval ; 
that is, if y(5) is a solution of (25) then 



-99- 


mU.-C J 

llyCcplI ^ ilyC5i)lle (105) 

where M < ~ can be determined for 1 < |C| < “ . There- 
fore, it is not necessary to continuously reduce (?) to 
canonical form in order to control the error growth. It 
is only necessary to reduce yj^(Q to canonical form 
several times over the interval of numerical integration. 
The length of the interval over Which a numerical integra- 
tion can safely be performed without reduction to canon- 
ical form obviously depends on the parameters cr, k, and y 
and the precision which is maintained in the calculations. 
To some extent experience is required, but a good first 
estimate of the length of an interval can be obtained by 
considering the different asymptotic rates of growth (see 
Section 4.1, step b in the numerical algorithm). 

The accuracy of the calculations, which were per- 
formed for the viscous problem formulated in Section 2, is 
discussed in Section 5 and Appendix B.2. 
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5. COMPUTATIONS AND CONCLUSIONS 

5.1 PRELIMINARY REMARKS AND ORGANIZATION OF THE 
COMPUTATIONS 


The viscous problem is solved for small c > 0 
when equations (32) and (33) are solved. The most diffi- 
cult numerical problem encountered is the determination 
of vectors at 5 = which specify the different asymp- 
totic solutions. There are several cases which must be 
considered separately. The individual cases can be clas- 
sified according to the character of the roots of the 
dispersion relation. 

CASE 1: The roots of the dispersion relation, 

and X^ > are real and distinct. In addition, 

— - 7 ^ 0 . 

Y 

Upon solving equation (32), constants e^, e^, e^, 

and e, are determined such that 
4 

(5) + 62^2(0 = inv 2(?) “ eg) 

- e^ TLSOL (g) . (106) 


If ^ 7 ^ ^2 ’ th,en 


+ ^25^2 (5) 



e _b^ 
3-0 




( 107 ) 
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2 1 

However, ^ becomes negligible as e 0 , 

for 0 < z < B < “ . For e ss 10~^^ the vector 

^ 2-^1 

®3— 0 1^ neglected below the transition region 

except for ^2 ” ~ ^ ‘ case a /y ~ k =0 is a 

special situation which is discussed in case 4. 

Thus, above the boundary layer and below the tran- 
sition layer, the solution of the viscous problem can be 
accurately approximated by a multiple of the inviscid 
solution with finite kinetic energy in an infinite column 
of fluid. 

If 0<A<z<B<“, then (107) implies that 
the solution with finite kinetic energy is approached uni- 
formly on this interval . The only calculations which must 
be performed for real and X^ are for X^ « X^ 

(case 3) and a^/y - k^ ~ 0 (case 4) . 

CASE 2: The roots of the dispersion relation (13) 

are complex. For this case the results of Yanowitch 
[2,3] and Lindzen [4] imply that for small 3 , the 
dimensionless, vertical wave number defined in (20), the 
reflection coefficient is large and for g large the 
reflection coefficient is small. Thus, it is expected 
that reflection is significant near the boundaries of the 
shaded regions in figure 1. 
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Aside from a scaling constant, the solution above 
the boundary layer and below the transition layer is 
approximately 


U(z)^ 


U^(z) 


U^(z) 

_W(z) 


W^(z) 

R 

W2(z) 


(108) 


where 


W. (z) 

X 


is a solution of the inviscid differential 


equation (12) for i = 1 or 2 and normalized so that 


U^(o) = 1 . In addition, the inviscid solutions are num- 



U, (z)' 


U„(z) 

bered such that 

Wj_(z) 

has upward energy flux and 

W2(z) 


has downward energy flux. The constant k is, of 
course, the reflection coefficient. 

If £ > 0 is specified, then is determined 

from the constant e^ in (32) . The scalar e^ might be 
considered the asymptotic reflection coefficient. The 


constant e^ is invariant (see theorem B2 in Appendix 
B.2) as e ^ 0 . However, k does change as e 0 . 

R 


If 5 = £n - 
e 


then 


I I is asymptotically invariant as e ^ 0 

R 

for fixed a and k since |e^| is invariant. 


(109a) 
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Arg ic ''2B6 is asymptotically invariant as 
R 


e 0 , for fixed o 


acoustic wave 
12 „ \ 


and k and the 


Y 


> 0 


C109b) 


Arg K +266 is asymptotically invariant as 
R 

e -> 0 for the gravity wave 
/ 2 „ \ 


Y 


- k" < 0 


(109c) 


Due to relations (109) , it is more useful to compute 
arg k„± 26<5 than arg k . 

R K 

It is reasonable to expect arg k to change as 
£ ^ 0 since the reflection takes place in the transition 
region. .he transition region is in the vicinity of 
z = %Vl and this point recedes to infinity as e +- 0 . 

Thus, the reflection coefficient can be expected to 
undergo a phase change as e 0 . 

In Appendix B (see B.2, relations (B74) through 

_ *7 

(B84)) it is shown that a fixed error tolerance (10 ) 

on the asymptotic solutions resulted in an approximate 
value of where 


n 

a 


2 


( 110 ) 
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and for |?| > tKe truncated asymptotic expansions 

-* 7 

have a minimum relative error less than or equal to 10 
For a fixed value of k > 0 , cr is bounded away from 
zero in both the acoustic and gravity wave regions (see 
figure 1, Section 1.2) as 3^0. If the initial vectors 
for the canonical numerical integration from to 

5 = i are computed by method II (see B.l), then it is 
obvious that the asymptotic initial vectors v/ith inviscid 
asymptotic behavior tend uniformly towards each other as 
3^0. If the vectors DC^(5), , INV^(g) , and 

TLSOL (g) remain linearly independent as 3 0 , then 

|K^(k,3)| ^ 1 (111) 

as 3^0 since |e^| in equation (32) tends to one as 
3-^0. In addition, if e > 0 is held fixed, then 

K^(k,3) -> -1 (112) 

as 3-^0 since e^ 1 . Relations (111) and (112) are 
an immediate consequence of (31) and the above assumption 
on the linear independence of DC^ (g) , DC^ (g) , INV ^ (g) and 
TLSOL (g) . 

CASE 3: The roots of the dispersion relation are 

equal, that is, ^ ^2 ~ 7 ~ ^ this case 

the solution of the viscous problem in the inviscid region 
is approximately 
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U(z)“ 

= a) 

U2(z)- 

+ B 

Ul(z) 

( 

W(z) 

I 

^ 2 ( 2 ) 


W^(z) 

) 


( 113 ) 


where 


U^Cz) 
w, (z) 


is an inviscid solution and, for 


i = 1 or 2 , is defined in (18) . Since the asymptotic 
connection relations (32) are invariant as e ^ 0 it can 
be shown that B °o as e 0 . However , 


B 


OUn |i 


CU4) 


as e -> 0 and hence 


U2(Z) 

W2(Z) 


cannot be neglected in the 


entire ir’^'iscid region. It should be noted that A 0 

as £ -> 0 since the kinematic boundary condition 

(W(0) = 1) requires AxB to remain finite as e 0 . 

CASE 4: The roots of the dispersion relation are 

2 

X = — and X„ = ^ or — - = 0 . In the 

1 Y 2 Y Y 

inviscid region the solution of the viscous problem is 
approximately given by 



/ 



'i 

^U(z)^ 

« C } 

-iJ.izy 

+ D 

rupz) ( 

/ 

W(z) 


^2 (Z) 


W^(z)Jl 


( 115 ) 
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whexe 


Ujz) 

W, (z) 


is an inviscid solution and, for 


i = 1 or 2 , is defined in (16) . The asymptotic connec- 
tion relations (32) determine the constant D and since 
the asymptotic connection relations are invariant it fol- 
lows that 


as e 
If Y 


|Dl = o(e ^ ^ 

0 . Hence D 0 as e 0 since 
1.4, then 


(116) 


J > ^2 


and 


« .714 (117a) 

« .286 . (117b) 


The kinematic boundary condition (W(o) = 1) requires 


W(z) = CW 2 (z) + CDW^(z) K 1 (118) 


near z = 0 as e 0 , but W^(z) is identically zero. 
Thus , 



(119) 


as £ -> 0 , that is, C-><» as e->0. A resonant 
situation develops as e 0 . For a finite value of e 
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like 10 there should be a noticeable resonant peak for 

0 ^ 2 

the scalar C in the vicinity of — - k = 0 . 

Y 

In order to describe the resonant peak, it is useful 

to calculate the modulus of C in relation (115) for 
2 2 

oVy - 0 , in addition to the obvious calculations 

required for o /y - k = 0 . For this purpose it is best 
to normalize the inviscid solutions so that solutions (16) 
are approached as a /y - k 0 . This is easily accom- 
plished; consider 


and 


- X . 

1 1 Y / 


U^(z) 


k[x + - 

\ 1 Y 

W^(z) 


1 

“U2(z)“ 


“ 1 ~ 

W2(Z) 


o^/y - k^ 



d(Y - 7)_ 


Xiz 


^ 2 ^ 


C120a) 


(120b) 


where and X^ are assumed real and 


^2 < 2 < \ 


( 121 ) 


for Y “ 1 • 4 . 

Not only is it important to organize the computa- 
tions but it is also necessary to consider the actual 
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computational devices which perform the calculations. The 
calculation of the DC solutions was accomplished on a 
UNIVAC 1108 in double -precision (about 18 significant 
figures) . Three hundred terms were generated from the 
recursion relation (AlO) in Appendix A and then summed at 
5 = i . Due to certain compiler difficulties it was not 
possible to perform all the computations on the 1108. The 
calculation of asymptotic initial vectors and the canon- 
ical numerical integration was carried out on an IBM 7094 
in single-precision (about eight significant figures) . 

The error check (step 5 of the algorithm in Section 4.1) 
indicated that approximately four significant figures are 
maintained in the calculations. 

5.2 NUMERICAL RESULTS 

This section contains the numerical results for 
the viscous problem formulated in Section 2. The calcu- 
lations were performed for e = 10~^^ , which is com- 
parable to a value in the earth’s atmosphere. The other 
dimensionless parameters, k and o , when equal to 
unity, correspond to a horizontal wavelength of about 
45 km and a frequency of 2.5 radians per minute, respec- 
tively . 

CASE 1: The roots of the dispersion relation are 

real and distinct. In Section 5.1 it was shown that the 
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only calculations required are for (case 3) and 

o^/y - r; 0 (case 4). 

CASE 2: The roots of the dispersion relation are 

complex or the inviscid solutions are wavelike in z . 

The reflection coefficient, k in (108), is computed for 
various values of a and k . Due to the remarks made in 
Section 1 about some recent research of Yanowitch, it will 
be worthwhile to compare and e or equivalently 

£n |k I and -ttB. To avoid confusion will denote the 

acoustic reflection coefficient and will denote the 

gravity reflection coefficient. In order to construct 
figures 2 through 7 (pp . 119-124) the values of k and 
Kr^ were computed for 3 = .01, .1, .2, .4, .8, and 1.6 . 
For k = .5 , K was computed for the additional 3 

R A 

values of .5, .6, .7, .9 and 1.0. 

If the reflection coefficients are considered as 
functions of the horizontal and vertical wave numbers, 
then figures 2 through 7 imply 

|KR^(k,3)l < kRG(k,3)| < e'""^ (122)^^ 

]T9 “ __ 

The calculations which were performed for 

k = .005 lead to the relation I ( • 005 , 3) I > e”^^ . 

Q RG 

However, |k (.005,3)1 - e was very small and about 

KG 

equal to the estimate of the error obtained from step 5 of 
the numerical algorithm. 



- 110 - 


In order to test relation (122) and determine any 
sensitivity to changes in y > some computations were per- 
formed for Y “ 4 . Surprisingly, the acoustic reflection 

■" IT B 

coefficient was in much better agreement with e than 

the gravity reflection coefficient for y = 4 . Thus, it 
appears that (122) is a quantitative rather than a quali- 
tative relation; nevertheless (122) is a useful summary 
of the calculations. 

In order to completely specify k and k it 
is also necessary to determine the argument of these com- 
plex quantities. Due to the remarks in Section 5.1 it is 
best to compute arg K^„(k,B)+235 and arg K„^(k,B)-236 to 
obtain a useful profile of the argument of the reflection 
coefficient as e -> 0 (tables I through VII). 

TABLE I. -ARGUMENT OF THE REFLECTION COEFFICIENTS 

FOR k = .005 


3 

arg k^^+236 

arg Kj^^-236 

0.01 

-3.12 

-2.99 

. 1 

-2.94 

-1.68 

. 2 

-2.67 

- . 515 

.4 

-1 . 71 

.983 

. 8 

1.71 

2.26 

1.6 

- . 567 

2.44 
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TABLE II. -ARGUMENT OF THE REFLECTION COEFFICIENTS 

FOR k = .05 



arg k^^+236 

arg K^^-2B6 

O 

O 

3.12 

-3.00 

. 1 

2.97 

-1.75 

. 2 

2.86 

- .641 

.4 

2.99 

. 788 

. 8 

-1.66 

2.03 

1.6 

-1.34 

2 . 20 


TABLE III. -ARGUMENT OF THE REFLECTION COEFFICIENTS 

FOR k = .25 


3 

arg K^^+236 

arg Kr^-236 

O 

O 

3.11 

-3.02 

. 1 

2.88 

-2.00 

. 2 

2.67 

-1.09 

.4 

2 . 53 

.099 

. 8 

-3.03 

1.14 

1.6 

1.17 

1 . 20 
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TABLE IV. - ARGUMENT OF THE REFLECTION COEFFICIENTS 

FOR k = .5 



arg K^^+2B6 

arg Kj^^-23'S 

.01 

3.10 

-3.11 

. 1 

2.76 

-2.79 

. 2 

2.54 

-2.41 

.4 

2.52 

- .308 

. 5 


.165 

.6 


.345 

. 7 


.438 

. 8 

-2.99 

.487 

.9 


.507 

1.0 


. 504 

1.6 

.547 

.127 


TABLE V. -ARGUMENT OF THE REFLECTION COEFFICIENTS 

FOR k = .75 


3 

arg k^^+236 

arg 

0.01 

3.10 

3.12 

. 1 

2.73 

2.89 

. 2 

2.54 

2.64 

.4 

2.68 

2.10 

. 8 

-2.61 

.804 

1.6 

.822 

- .664 
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TABLE VI. - ARGUMENT OF THE REFLECTION COEFFICIENTS 

FOR k = 1.0 


B 

arg k^^+236 

arg k^^-236 

0.01 

3.11 

3 . 10 

. 1 

2 . 83 

2 .75 

. 2 

2 . 74 

2.38 

.4 

3.06 

1 . 74 

. 8 

-1.97 

.665 

1.6 

1 . 55 

-1 . 24 


TABLE VII. -ARGUMENT OF THE REFLECTION COEFFICIENTS 

FOR k = 1.5 


3 

arg k^^+2B6 

arg 

0 . 01 

3.13 

3.08 

. 1 

3.07 

2 . 56 

. 2 

-3.07 

2.00 

.4 

-2.32 

1.03 

. 8 

- .441 

- .409 

1.6 

-2.67 

-2.59 
















-114- 


The figures 2 through 7 and tables I through VII 

can be used to construct a multiple of the solution to 

the viscous problem above the boundary layer and below 

the transition layer. Thus, it is desirable to obtain 

some estimate of the lower boundary of the transition 

layer and the upper boundary of the boundary layer. 

In Section 3.2, relation (30) implies that the 

boundary layer solution decreases by a factor ^ as z 

/2 yg - 1 1 

is increased by e = 10 a linear combi- 

nation of inviscid solutions agrees with the solution to 
the viscous problem to four significant figures if 

"b.i. ^ i 


"b . £ 


6 X 10' 

/a 


(123) 


The upper boundary of the boundary layer varies from 
several meters for the gravity wave (k = .005, 3 = 1.6) 
to less than one-half meter (k = 1.5, 3 = 1.6) for the 
acoustic wave. 

In order to obtain a more detailed description of 
the solution to the viscous problem, some computations 
were performed which illustrate how the transition region 
joins the inviscid and viscous regions. It is useful for 
this purpose to consider plots of | U (5) | x/ | ^ | (ordinate) 


versus £n /\K\ (abscissa) and | W(^) | x/ | ^ | versus £n /JX' 



-115- 


The ordinates |U (?) 1 x/lCj ] W(^) ] /| 5 [ proportional 

to the square roots o£ the kinetic energies associated 
with the horizontal and vertical velocity components, 
respectively. The abscissa, £n /jTT, is proportional to 
z/2 and hence a one-unit change in £n / ] C 1 corresponds to 
two scale heights (see figures 8 through 14) . It should 
be noted that convenient multiples of the viscous solu- 
tions are plotted in figures 8 through 14. 

In figures 8 through 14 the inviscid region is 
easily identified. It appears that the inviscid region 
begins at about z = 20 (140 km) and extends downward 

from this height, that is, approximately 20 

scale heights. The oscillation about the dotted line is 
caused by alternate constructive and destructive inter- 
ference of the inviscid solutions with upward and down- 
ward energy propagation. The dotted line, in figures 8 
through 14, corresponds to a multiple of the inviscid 
solution which has upward energy propagation and hence 
satisfies the radiation condition. The amplitude of the 
oscillation about this dotted line is an indication of 
how poor or good an approximation the radiation condition 
is. For 8 = 1.6 the radiation condition appears to be 
quite good (in the inviscid region); if 8 = .4 it is 
unsatisfactory. Relation (122) implies that the 
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radiation condition introduces an error o£ about 4 per- 
cent for B = 1 , and for larger values of 3 a 
smaller error is expected. Thus, if the vertical wave- 
length is less than 45 km (3 > 1) the radiation condition 
is substantially correct. 

CASE 3: The roots of the dispersion relation are 

1 1 

equal or ci = and ~ ^2 ~ constants A 

and B in relation (113) are listed in tables VIII and 
IX. 


TABLE VIII. -THE CONSTANTS A AND B , DEFINED IN 
RELATION (113), FOR o^/y - k^ < 0 AND y = 1-4 . 


k 

o 

A 

B 

■ 

arg B 

0.005 

0 . 00535 

-11. 5-. 935 i 

-25.4+1.57 i 

25.5 

3.08 

. 05 

.0534 

-1.07-.078 i 

-27.3+1.57 i 

27.3 

3.08 

.25 

. 260 

-.170-. 013 i 

-27 .7+1.69 i 

27.8 

3.08 

. 5 

.447 

-.044-. 0036 i 

-28.5+1.87 i 

28 . 5 

3.08 

.75 

. 505 

- . 018- . 0013 i 

-28.7+1.66 i 

28 . 7 

3.08 

1 . 0 

.520 

-.012-. 00083 i 

-28.2+1.58 i 

28 . 2 

3.09 

1.5 

. 529 

- .0073- .00054 i 

-27.0+1.57 i 

27 .0 

3.08 
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TABLE IX. - THE CONSTANTS A AND B , DEFINED IN 
RELATION (113), for a^/y - > 0 AND y = lA . 


k 

o 

A 

B 

■ 

arg B 

0.005 

0.592 

0 .00015- . 0000084 i 

-34.0-1.57 i 

34.0 

-3.10 

. 05 

. 592 

. 00154- . 000088 i 

-33.6-1.59 i 

33.6 

-3 . 09 

. 25 

. 608 

,0099-. 00084 i 

-32.3-2.27 i 

32,4 

-3.07 

. 5 

.707 

. 043- . 0083 i 

-28.1-4. 31 i 

28.4 

-2.99 

.75 

.939 

.120- . 015 i 

-25.1-2.46 i 

25.2 

-3.04 

1.0 

1.22 

.199-. 019 i 

-24.4-1.80 i 

24.4 

-3.07 

1.5 

1 .79 

.357-. 032 i 

-23.4-1.59 i 

23.4 

-3 . 07 


CASE 4: The roots of the dispersion relation are 

X = i and X = or o^/y - = 0 . In addi- 

1 Y 2 Y 

tion, computations were performed for 

= (1 ± .01)/7 k . (124) 

The results of the computations are summarized in 
figure 15. 

There is a noticeable resonant peak for 

2 2 

o / Y - k = 0 . Varying a by 1 percent results in a 
decrease of |C| by approximately a factor of 10 (refer 
to relation (115) and figure 15) , 
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The boundary layer solution is not capable of 
reducing U(z) in (115) to zero in 0 < z < z^^ ^ with 
negligible effect on W(z). The value of |C| computed for 
e = 10 depends on the peculiar properties of the 

boundary layer solution. Other effects in the earth's 
atmosphere more important than viscosity are being 
neglected near z = 0 , Therefore, the results for this 
case are of dubious value . 
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0 2 4 6 8 10 12 


)ln S 

Figure 8. — Viscous gravity waves for k = .05 and B = 

s =/R| • 


.4, 1.6 
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Figure 9. — Viscous gravity wave s fo r k = .5 and g = 

s = /rtf . 


.4, 1.6 
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Figure 10. — Viscous gravity waves fo r k = 1.5 and g = .4, 1.6 

S = /U1 . 
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Figure 13. — Viscous acoustic wav es f or k = .5 and 3 = .4. 1.6 

s = /UI • 
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Figure 14. — Viscous acoustic waves for 

s = yul 


k = 1. 5 and 3 


.4, 1.6 




An k 


Figure 15. — The resonant Lamb wave and 1 percent variation in 

the frequency for e = 10”^^ and y = 1-4 . The scalar C 
is defined in relation (115) . 
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5.3 CONCLUSIONS 


For the case of real distinct roots of the disper- 

2 2 

sion relation (13) and a /y - k 0 , the solution of 
the viscous problem approaches a multiple of the inviscid 
solution with finite kinetic energy in a column of fluid 
of finite cross section as e -5- 0 or p 0 . Moreover, 
the convergence is uniform on any interval 
0<A<z<B<“ as y-^0 (see Section 5.1, case 1). 
This conclusion is reasonable and in complete agreement 
with the results obtained by Yanowitch [2,3], If 
a /y - k =0 a resonant situation develops and the 
limiting case, as e ->■ 0 , results in a viscous solution 
with a horizontal velocity amplitude which is 


which is 0 


llY 
, e 


viscosity therefore 
in Sections 5.1 and 
By far the mo 
roots of the dispers 
tions which are wave 
Yanowitch [2,3] and 
expect that 


near z = 0 . A nonzero value of 

limits the resonant peak (see case 
5.2) . 

st interesting case concerns compl 
ion relation, that is, inviscid so 
like in z . Due to the results o 
Lindzen [4] it seems reasonable to 


4 


ex 


lu- 

f 


-IT 


R 


(125) 
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However , figures 2 through 7 indicate that (125) 
is not satisfied for all horizontal wavelengths. Thus, 
the modulus of the reflection coefficient depends on the 
model selected; for example, the inclusion of compressi- 
bility appears to be quite significant. 

A careful examination of figures 2 through 7 indi- 
cates that f(k,B), where 

f(k,3) = |tt3 + iin |k (k,3)|| , (126) 

K. 

has a maximum near k = 3 . If more general problems are 
considered, it is probably true that the mathematical 
model will play a major role in determining at 

least for small values of 3 . 

Several qualitative statements can be made for the 
case of solutions which are wavelike in z . 

a. Reflection is negligible for 3 large, that 
is, the radiation condition is substantially 
correct (refer to discussion of case 2 in 
Sections 5.1 and 5.2). 

b. Reflection is important for 3 small and the 
reflection coefficient tends to the limiting 
value -1 as 3-^0 for fixed e > 0 and 

k > 0 (refer to discussion in Section 5.1 
and relation (111) ) . 
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c. The modulus o£ the gravity reflection coeffi- 

“ 7T 6 

cient more nearly equals e for all values 

of k than does the corresponding acoustic 
reflection coefficient for y = 1-4 (see 
figures 2 through 7) . 


The problem formulated in Section 2 is limited to 
small oscillations of the lower boundary z = 0 , An 
estimate of the validity of the linearization can be 
obtained by determining the maximum value of 

is a solution of the vis- 



U(z)l 

, where 

'U(z)' 


W(z) J 


W(z) 


cous problem. For the cases which were computed, the 


maximum of 


U(z) 

W(z) 


satisfied 



U(z)^ 



W(z) 



4 X 10' 




( 127 ) 


20 


where the boundary layer is given by 0 < z < z „ and 


"b . il . 


,-11 


£ is assumed to have the value 10 

The kinematic boundary condition, W(o) = 1 , is 

H 

equivalent to requiring an amplitude of — of the 


^Estimate (127) is correct for k = .05, 0.5, 1.5 , 
and B = .4 and 3=1.6 for both the acoustic and 
gravity waves. For real distinct roots of the dispersion 
relation, a considerably smaller bound should result 
since the inviscid solution with finite kinetic energy in 
an infinite column of fluid grows less rapidly than 
z/ 2 

e as z increases. 
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oscillation o£ the lower boundary. The araplitude o£ the 
oscillation scales the solution of the viscous problem. 
Hence, an oscillation of amplitude 10 ^ H or 0.7 cm would 
be more reasonable, that is, consistent with the lineari- 
zation. Of course the resonant case a /y - k =0 
would require a much smaller amplitude in order for the 
linearization to be valid since |UCz)j is very large. 
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APPENDIX A 

THE REGULAR SINGULARITY 5=0 

In this appendix the regular singularity 5=0 
is investigated. Two linearly independent solutions of 
the viscous differential equation (25) are developed 
about 5=0, which satisfy the DC (11) . In addition, 
it is shown that there exist two other solutions of 
equation (25) which violate the DC, 

A.l THE DC SOLUTIONS 

Assume that a fundamental matrix solution of 
equation (25) can be expressed in the form 

S 

m 

where 

S are constant square matrices for (A2a) 

m = 0 , 1 , 2 , • • • 

J is the Jordan canonical form of R (A2b) 

(defined in relations (24)) 


$(?) 


E 

m=0 


-m 


(Al) 


2 1 

In all subsequent contexts the functions £n5 will 
be defined to be the principal branch. 
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Sq is a constant nonsingular matrix such (A2c) 

that SqJ = RSq 

principal branch (A2d) 

of £n^ 


In order to construct the eigenvalues and 

eigenvectors of R will be determined. It is easily 
shown that the eigenvalues of R are +k and -Ic and each 
eigenvalue is of multiplicity two. Only one-parameter 
solutions are obtained upon solving 



(R - kl)e^ = 0 

CA3a) 

and 

(R + kl)e.3 = 0 , 

(A3b) 


where 



let a = 1 (A4a) 


and 



let b = 1 . (A4b) 


Thus, the Jordan canonical form of R is 


J 


k 1 0 0 

0 k 0 0 

0 0 -k 1 

0 0 0 -k 


(A5) 
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In addition to the vectors e^ and e^ the generalized 
eigenvectors e^ and e^ (solutions o£ (A6)) must also 
be determined. 


(R - kDe^ = e^ 


(A6a) 


and 


(R + kl)e^ = £3 


(A6b) 


If the vectors e^ and are normalized by setting 

the first components equal to one, then 


and 


1 

k + 1 

k - 7 
k 

k - 6 


1 



(A7a) 


(A7b) 


k + 6 
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The vectors e ^ define the matrix 

that is , 


1 

k 

1 

k 


1 

k + 1 

k - 7 
k 

k - 6 


1 

■k 

■1 

k 


1 

1 - k 
k + 7 


k 

k + 6 


(A8) 


If CAl) is substituted into (25) and is 

defined by (A8) , then the recursion relation 


+ 1)1 + J] + KS^[inI + J] 


RS + DS (A9) 
m+1 m 


is obtained. If denotes the ith column of S , 

— m m ’ 

then (A9) can be written as 


(m + 1 + Jl ,)s + 6 ,s + 

^ m+1 1 -in+l 


K 


(m + £,)s^^^ + 6,s 

3 -m 3-in 


(j-l) 


Rs^^J + Ds^^^ 

-^n+l — m 


(AlO) 


where 


0 , 6 , 


1 (Alla) 


and 


k , 1 . 


-k . (Allb) 


The recursion relation (AlO) can be placed in the form 


s (j) 
— in+ 1 




m -in 


D^i) /s(:i-i 

3 m \— m+1 


+ Ks 


(j-l) 


(A12) 
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where D ^ . + lU + 1) I - R] ^ (A13a) 

m 3 

and = [a, + m + 1)1 - R]"^[D - (m + il.)K] . 

(A13b) 

Due to the DC CH) the vector solutions corre- 
sponding to j = 3 and 4 will be discarded. Therefore, 
it is only necessary to develop the recursion relation 
(A12) in detail for j = 1 and 2 . 

For j = 1 or 2 , = k . Thus, consider 
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and 


[D - Cm + k)K] 


0 


0 0 0 



0 


k 

Y 


k 


0 


0 0 0 


-3(y - l)k 
4y 


3k 3a^ 3g 

"T" "TY" T" 


(A14c) 

where g = k + m + 1 (A14d) 

It is clear that the recursion relation (A12) is 
well defined as long as the matrix [ (£ ^+m+l) I-R] is 
nonsingular for m = 0,1,2,*** . Recall that [-kl-R] 
and [kl-R] are singular matrices since k and -k are 
the eigenvalues of R . Moreover, k is the only 
nonnegative eigenvalue of R and, hence, [(k+m+l)I-R] 
can never be singular for m = 0,1,2,* •• . Thus, the 
recursion relation (A12) is always well defined for 
j = 1 or 2 . 

Quite obviously, if k is an integral multiple 
of one-half, then the recursion relation (A12) breaks 
down for j = 3 or 4 since [ (-k+m+1) I -R] is singular 


2 2 

The recursion relation (A12) is well defined if 
the vectors , for j = 1,2,3 and 4 uniquely 

specify the vectors ^ for j = 1,2,3 and 4 and all 

positive integral subscripts m . 
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when m = 2k - 1 . This difficulty can be overcome by 
introducing a shearing-type transformation [8, Chapter 
4, Section 4] of the differential equation (25). This 
tactic is unnecessary since the solutions corresponding 
to j = 3 or 4 are to be discarded because they violate 
the DC. The approach adopted for j = 3 or 4 is one 
that is used over and over again in asymptotics and is 
discussed in Section A. 2. 

Since the recursion relation (A12) is well defined 
for j = 1 or 2 the first two columns of 4(0, defined 
in (Al) , are solutions of equation (25) . Denote the 
first two columns of 4(5) by DC (5) and DC^ (5) . Thus, 

00 

DCi(C) = 2 (A15a) 

m=u 

oo 

and + (£n 5)^^(0 . (A15b) 

As 5 0 , DC (5) and DC^ (^) exhibit the scalar 

growths 5^ and (£n 5)5^ o^- e'”^^ and ze~^^ as z -> «> . 

Hence, DC^(5) and satisfy the DC (11). The 

expansions (A15) converge for 0 < |5l < 4/3 since the 
nearest nonzero singularity occurs at 5 = 4/3 . The 
solutions DC^(5) and DC^ (5) can be analytically continued 
to the whole ray | arg 5 = j since the differential 
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equation (25) has no singularities for 0 < |5| < ~ and 
arg ? = . 

It should be noted that in addition to the decay 
and (iln 5)?^ the solutions DC^ (C) and DC^ (g) are not 
wavelike in z for sufficiently large z or equiva- 
lently for sufficiently small ^ ; that is, the 
variation in the argument of the components of DC (g) 
and DC^ (g) is bounded in the vicinity of g = 0 . This 
is a consequence of 

arg(g^) = constant (A16a) 

and arg(£n g) -> constant (A16b) 


as 5 0 for arg 


1 

arg i: — constant 


im=0 


(A17) 


CO 

as g -> 0 ; that is , the components of E 

m=0 

to constants as g -> 0 . Thus, in the viscous region 
the wave motion is damped and not wavelike in z . 


A. 2 THE SOLUTIONS WHICH VIOLATE THE DC 


In this section the vector solutions of (25) which 
violate the DC are investigated. If k is not an inte- 
gral multiple of one-half, then the recursion relation 



-145- 


(A12) is well defined for j = 3 and 4 . The third and 
fourth columns of 4(5) defined in (Al) , exhibit the 
scalar growths 5 ^ and (£n 5)5 ^ as 5 0 or e^^ 

and ze^"^ as z . Hence the third and fourth columns 

of 4(5) violate the DC (11). If k is an integral 
multiple of one-half, then the third and fourth columns 

of 4(5) cannot be constructed from (Al) since the 
recursion relation (A12) is not well defined. However, 
it can still be shown that there exist two solutions, of 
the differential equation (25) , which exhibit the scalar 
growths 5~^ and (£n 5)5~^ as 5 0 . Thus, for all 

values of k there exist precisely two linearly inde- 
pendent solutions which satisfy the DC (11) and two 
other solutions which violate the DC. 

Consider 

4(5) = Sq5‘^ , (A18) 

where 

5)e^)5’',e35“’",( 

(A19) 

and , £ 2 ) £3 3 .nd e^ are the column vectors which 
define (see (A 8 )). It is easily verified that 

g|i(C) = (1) ho . (A20) 


e^ + (£n 5 ) 63)5 
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Suppose 4(5) is a fundamental solution of 
equation (25) , then 

^$(5) = (|+0(l)j4(5) . (A21) 

The symbol 0(1) is used in (A21) to denote that the 
remaining description of the coefficient matrix is 
bounded as 5 0 . 

/\ 

Suppose the ]th column of 4(5) is denoted by 
$^( 5 ) and the ]th column of 4(5) is denoted by ^^(5)- 
If the first two columns of 4(5) are chosen to be the 
convergent expansions (A15) , then it is easily verified 
that 

Ili.(5) - ij(C)ll = 0(1) X Ili.(5)|l as 5^0 , 

(A 22 )^^ 

for j = 1 or 2 . If it can also be shown that there 
exist solutions ^ 3 ( 5 ) and ^^(5) of equation (25) such 
that relation (A22) is valid for j = 3 and 4 , then 
it immediately follows that ^ 3 ( 5 ) and ^^( 5 ) violate the 
DC (11). Thus, 4(5) would consist of two column vectors 
which satisfy the DC and two solutions of (25) which 
violate the DC. The remainder of this section is devoted 

^^The maximum norm is used throughout this paper. 
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to establishing the existence of solutions ^ 3(0 and 
^^(5) which satisfy (A22) . 

The method of attack is to treat the differential 
equation (A21) as though the term 0 (l)x$ (^) , on the 
right-hand side of (A21) , is an inhomogenous term of 
the differential equation. The homogenous part of the 
right side of (A21) is ll-j , that is, the homogenous 

differential equation is (A20) . If the term 0 (l)xf( 5 ) 
were a known function, then a closed-form solution of 
the differential equation (A21) could be obtained by 
the method of variation of parameters. Since 0 (l)xf(^) 
is not known, an indefinite integral equation is obtained. 

Attempt to determine a vector solution j of 

equation (A21) , of the form 

l(^) = S(?)c(C) . (A23) 

If equation (A23) is substituted into the differential 
equation CA21) , then 

i(?) = J^$CC)i“^(S)0(l)^(S) dS CA24) 

is obtained, where 0(1) is defined in (A21) . 

If 

$^(S) + $^(S) = i“^(S) 


(A25) 
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and ij(5) 


is a continuous solution of 


4). (5) = i.(C) + (s)o(i)0 (s) ds 

3 3 O 3 

+ J^$(a^,(s)o(i)(f) rs) ds , 

*a ^ ~3 

then is a differentiable solution of 

equation (A21) . Let 


$^(S) 


r ~ 

0 

0 

0 

0 


0 0 
0 0 
0 

0 0 


0 

0 

- (£n S)S^ 



and 


S"^ - (£n S)S"^ 0 

0 S“^ 0 

0 0 0 

0 0 0 


0 

0 

0 

0 


-1 

0 


Now consider the Picard iteration 

= 0 

and 

= $.(^) + (S)O(l)*^^ (S) dS 

~J — 3 0 3 

+ /^$(?)$,(S)0(l)i_|'^^ (S) dS , 
a ^ 3 


(A26) 


(A2 7a) 


(A27b) 


(A28a) 


(A28b) 


(A28b) 
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where 0 < |^| < |a| and a, and S lie on the ray 

arg S = -j and j = 3 or 4 . 

LEMMA A1 : The Picard iteration, defined in (A28) , 

converges for |a| > 0 sufficiently small. The con- 
vergence is uniform for each compact subset of the line 
segment 0 < < |a| and arg ? = addition, 

if |a| > 0 is sufficiently small, then the iterants 
satisfy the inequality 


where 

(116311 if 3=3 
m.(?) = { ^ 

j I le^ I I + Un 5 1 I |e^ I I if j = 4 


(A29a) 


(A29b) 


PROOF: (induction) 

Thus, for n = 0 the relation (A29) is satisfied. The 
lemma will now be established for all values of n by 
induction. Assume relation (A29) holds for some 


n 


1 > -1 


and consider 
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Mi] 


(n+l) 


CO - 4) 


CO I I < ( 

3 0 -*- \ 


£n 


+ 1 


I <}> 


(n) 


- 4, 


(n-1) 


CS) I I X |ds 


0O( 


£n 5 


+ 1 


I i t 


- i- 


(n-1) 


CS) 1 1 X |ds 


where M and M are chosen so that 

X ^ 


1 |$CS)f^CS)0CD I 1 < 


5,n 


+ 1 


for 0 < |S| < |?| < [a] for some |a| >0 . 


Similarly 


$(0$2(S)ocD 1 I < M. 


in I 


+ 1 ’ 


for 0 < 1^1 < |S| < |a| . The inductive hypothesis 
be applied to the integrand in (A30) . Thus 


ll^y+i)(E) . jy>(5)ii < / 

J J o \ 


in I 


+ 1 


k /^\ n 


X m 


(0|dS| + fh\ { in I 


CS) 


"^CS) 

(A30) 


can 



CA31) 
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The first integral on the right side of 

o 

(A31) has a singular integrand (logarithmic growth 
near S = 0) . However, a logarithmic singularity is 
very weak and any power of £n S is integrable on a finite 
portion of the ray |arg S = -jj passing through S = 0 . 
Thus, the integral from zero to ^ in (A31) satisfies 

P = 0(1) X 1^1 X ^ (A32) 

as I a I •> 0 . 

The integral from a to C in (A31) can be 
bounded by 




n-l 


M. 


-k 


P\ 


In 


+ 1 


k 


m. (S)|dS 


(A33) 


However , 
£n ■ 


tends to zero faster than any power of 


tends to zero. Thus, 


tends to infinity as 
the integrand on the right side of (A33) is uniformly 
bounded for 0< |C|£|s|<|a( . Therefore, 




0(1) X 5 


n-l 


(A34) 


as 


a I ->■ 0 . 
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If |a| > 0 is chosen sufficiently small, then 
relations (A32) and (A34) imply 

p + < (I) X |5|-^ X m (^) . 

o a \ ^ / J 

NOTE: The value of the constant a can be 

chosen independent of the superscript n since a must 
only be chosen sufficiently small to satisfy 



for 0 < 1^1 £ |a| and clearly this can be done since 

E E 

the constants and and functions £n ^ ^ 

and nijC?) do not depend on the superscript n . 

Thus, relation (A29) is established by induction, 
but (A29) implies 

M+K+1 , , 

X lEr’^m. CE) 


if ^(5) 
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since 

+ ({)-♦<“> (5)1 1 < I (C) - col I 

- llif"’^>(0 - 4“"’^-^>(E)|| 

t ... + I (5) 

-l‘"’(E)ll . 


Hence, the Picard iteration converges uniformly in ^ on 
each compact subset of 0 < |c| < |a| and arg ^ . 

Q.E.D. 

THEOREM A1 : The Picard iterants and 

^ (5) j defined in relations (A28) , converge uniformly 
to the solutions and <|).(?) of the differential 

— j -i-4 

equation (25) , respectively, and 


I li, (5) - .$..(?)! I 


0(1) X 




(A3 5) 


as 5 0 and arg ^ = -j . 


PROOF: Since the Picard iterants ^^’^^ (^) converge 

uniformly for each compact subset of 0 < |^| <. |a| 


and arg 5 


TT 

'Z ’ 


thus the limit function ^^ (?) , where 


ij(?) 


lim (n)( ) 
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is a solution o£ the differential equation (25) for 
0 < |Cl <• |a| . This follows because (C) is a con- 
tinuous solution of the integral equation (A26) and, 
hence, is a differentiable solution of the differential 
equation (25). In addition, the estimate (A29) implies 

|lij(OM < 2Ur’"m,(0 , 

where m^ (?) is defined in lemma A1 . Thus, 

lllj(53 -IjCOll < £n I + 1 j I ^Isr^m. (S) |dS| 

• |dS[ . (A36) 

The integrals on the right side of (A36) must be 
estimated separately as ? ->- 0 . 

< Ur’" /^2M^(^|iLn || H- ljm.(S)ldS| = o(l) x 1 5 | 

since any power of £n ^ is integrable near ^ = 0 . 

The integration from a to ? on the right side of 
(A36) is not handled so easily. For |?| >0 suffi- 
ciently small the following inequality is satisfied: 

0 < I? I < ^/TTT < la| . 
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Consider on the right side o£ (A36) ; then 

a 



dSl 


(A37) 


The integrand on the right side o£ (A37) is bounded as 

k , 


5 0 since 


£n 


+ 1) is bounded. The £actor in 


£ront o£ the integral on the right side o£ CA37) is 
obviously o(l)x|^| ^m.(5) as |Cl 0 . 


All that remains to be shown is that 



(A38) 
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but the integrand on the right side of (A38) 


bounded for 0 < 


S - ^ 


The integrand is 
over an interval which shrinks to zero as ^ 
Thus , 



0 ( 1 ) X |c|-^m.(5) as |5| 


is uniformly- 
integrated 

0 . 

^ 0 . 


Q .E .D . 
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APPENDIX B 

THE IRREGULAR SINGULARITY ^ 

In this appendix the formal asymptotic expansions 
about 5 = CO are developed and investigated. Two 
methods are considered. One method of obtaining asymp- 
totic solutions, referred to as method I, consists of 
transforming the differential equation (25) until a guess 
can be made concerning the form of a fundamental set of 
asymptotic expansions. The other method, called method 
II, avoids the many transformations required in method I. 
However, method II only determines the asymptotic expan- 
sions with algebraic growth in ^ , that' is, only the 
asymptotic expansions with lead terms which correspond to 
the inviscid solutions are determined. Multiples of the 
boundary layer and transition layer solutions can be 
determined by numerical integration. Thus, for large 
values of 5 a fundamental set of approximate solutions 
of equation (25) can be obtained by method II. 

The proof of the existence of actual solutions of 
the differential equation (25) which are asymptotic to 
the formal expansions about the irregular singularity, 

C = °° , follows a rather standard format [8, Chapter 5]. 
For this reason no attempt will be made to establish 
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the asymptotic nature of the formal expansions, but this 
property will be used freely. 

B.l FORIdAL ASYMPTOTIC EXPANSIONS 


In this section the formal expansions of the dif- 
ferential equation (25) about the irregular singularity 
^ = oo are developed for four distinct cases: 

CASE 1: The roots of the dispersion relation (13) 

are distinct. In addition, 2(X^ - ^ integer and 

a^/y - 0 . 

CASE 2: The roots of the dispersion relation are 

equal or ~ ^2 ~ 7 ‘ 

CASE 3: The roots of the dispersion relation 

equal an integral multiple of one -half. 

CASE 4^^: CT^/y - k^ = 0 . 

CASE 1: 2(X^ - X^) f integer and o^/y - k^ / 0 . 

Consider the differential equation (25) which can 
be written in the form 





> 


where A ^ ^ converges for 1 5 | > 4- and 

O n ^ 


(Bl) 


^"^Method I is used for cases 1, 2, and 3; method II 
is used for case 4. 
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A 


Y 


- 


0 

0 


- ki 

0 
0 


(B2a) 


A_ 


4k^ 

3 

0 

Y - 1 
Y 


0 

O' 

k 


0 

0 


Y 


(B2b) 


and for ' n > 2 


A = 
n 


0 
0 
0 

4\^"^ Y - 1 


Y 


0 

0 

0 

n- 2 


0 

0 

0 

r'lif ^ 


n-l 


(B2c) 


2 

Since A^ = 0 , equation (Bl) falls in the class 
of differential equations v;ith a nilpotent lead coeffi- 
cient matrix. The formal asymptotic solutions are not 
trivially determined for this class of differential equa- 
tions. One approach is to transform (Bl) so that the 
transformed equation has a form for which an immediate 
guess can be made concerning the asymptotic solutions 
[10], This approach (method I) will now be developed. 



The first transformation to be considered is a simi- 


larity transformation. Let 

iCC) = e£(c) 


(B3a) 


Inhere 



(B3b) 


CB3c) 


The vector fiX) then satisfies the differential equation 
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where 



(BSb) 


and for n > 2 
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0 
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0 


0 
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0 


0 

A * 
n 

-&r 

Yf) 


ar 


1 

1 

(?n j ^ - WSr • »■) (if 

Yf) 


an 

- t) 






(B5c) 

and 

the scalar a is defined in (13) . 





Now consider the shearing transformation 



2_C?) = diagonal |l, 1, 

t) 

y (t) , 


CB6a) 


H 

ar g ^ 




where x = / 1 5 [ e 

2 



(B6b) 


1 

~1 0 

0 o' 



and diagonal^l, 1 , 1 , ^ j 

0 1 
0 0 

0 0 
1 0 

j 

CB6c) 



0 0 

» i 




The vector satisfies tke differential equation 


dy (x) 
dx 



CB7a) 





CB8b) 
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CB8c) 

and for n > 2 



(B8e) 
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Now consider anotJier similarity transf ormat ion , 


Let 


yCx) = Fy_(x) 


(B9a) 


where 


F = 


1 


Y 


Y - 1 
0 
0 


a 


1 

0 

■ya 

0 


0 0 


0 0 


a 


o 


1 1 


(B9b) 


and 




aa 


- Y - 1 
2 

Y a 


Y - 1 
2 

Y a 


— /y g (y - 1) 


2y 

/y g(Y - 1) 

o 2 

2y 


0 0 


0 0 


o 


2 /y 


2/7 


(B9c) 


then yC^) satisfies the differential equation 
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where 


Thus , 


dy (t) 
di 



(BlOa) 


1 = F . (BlOb) 

n n 


> 


(Blla) 

/y 

’ /7/ 



2U-I) 

Y 


-2/y(y- 1) 2/7tY-I) ■ 


2 2ya 

T ■ TwT 


2(y-D 

Y 


z-^Cy-1) 

Y^oci 


-2/y(y-I) 


^ 0 ^ 4k 


Y ' 1 




■ - ity 


b (y - 1) 




2o + 2 3 

Y ■ 2 


-(3,1) ELEMENT 


-(3,2) ELEMENT 


2o‘ ♦ 2 3 

— y I 


CBllb) 
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0 


-2H 


-f-T " ^ 

yk -a^ Y^a 


G + ^ Ck^ - o^)a G - yk^a ^5^ » ^ H °tT-l) + 


a h a^*l\ 

Y- 


+(3.1) ELEMENT +(3,2) ELEMENT 


k7 k /? °a ■ „ /f g(T-l) 

O .2 2 2 

Yk -o Y 


o l2 o^ + l\ 


and for n > 2 
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1 


2n-l 


0 
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/4\""^(o^-k^)Ya 

-TFTj^ 



[kV 


G] 



-(3,1) ELEMENT 


-(3,2) ELEMENT 



+ U,3) ELEMENT 


+(2,3) ELEMENT 


•(4,3) ELEMENT 


'(3,3) ELEMENT 


CBllc) 
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0 
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(l) tG - k^Y“l 
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(Bile) 

where 






G 

= 4 Y 

3 

■ 1 1 2 . 2 

k + (a 

Y ^ 


(Bllf) 


H 

= 1 + 

1 


(Bllg) 


_2 ,2 



Sibuya-type transformations are now considered. 
The purpose of these transformations is to reduce the 
off-diagonal elements of the coefficient matrix to zero. 
This procedure will be modified. Only the off-diagonal 
elements of the matrix associated with — will be 

T 

reduced to zero. This is a simpler task and ensures the 
convergence of the transformed coefficient matrix [10, 
Chapter IV] . 

Let 

y (t) = 


P(t)X(t) 


(B12a) 
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wher e 


P (t) 



(B12b) 


Pq = I (identity matrix) , (B12c) 


and 


/ y-1 \ / y- 1 \i 

(ya^aj (ya = J: 



y^c= . l' 

, 4k^ 
a + -3- - 

, . - ip 

o^a 

Q 

^s> 

K) 

0 

0 


2(Y-1) 2 

0 ^ (y-l) 


2 


+ (3,1) ELEMENT 0 0 


(B12d) 


then X(t) satisfies the differential equation 


dX(T) 

dr 



where 


C 


0 


2a la \ 

? > I 

/y /y/ 


(B13a) 


diagonal (0 , 0 


(B13b) 
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-i 

'y 


2ya 

Y-1 


2Cy-1) 

Y 

2(y-1) 

Y 


1 

2 


2a^ + 2 
Y 


3 

2 


2a^ + 2 


3 

2 


CB13c) 




(B13d) 


and for n > 3 


C = A + A P - P C , 
n n n-1 1 1 n-1 


(B13e) 


One additional transformation, a combination of 
a similarity and a Sibuya transformation, reduces to 

diagonal form and leaves unchanged. Consider 

X(t) = P(t)Z(t) , (B14a) 


where 


P(t) 



CB14b) 



"Td i 
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^lY 


Y - 1 
0 
0 


1 


V 

Y - 1 
0 
0 


0 0 

0 0 

1 0 
0 1 


(B14c) 




1 - X^Y 


Y (>^- 




Y) 


Y - 1 




Y - 1 


yC^^ - x^) 


0 

0 


0 0 

0 0 

1 0 
0 1 


(B14d) 


and 


0 0 


0 0 


0 0 


0 0 1 
4a\ Y 2, 


3 

Ta \2 


2a^ + 2 
Y 


(B14e) 


RECALL: and are the roots of the dispersion 

relation (13) , 
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Thus , 


dZ(x) 

dr 


Vn=0 / 


where 


Bq = diagonal/o, 0, — , , 

\ /y /y / 


(B15) 


CB16a) 


= diagonal ^-2X^, j, 


1 1 
2 ' ’ 


CB16b) 


^ = ^0 ^ " Pi> > 


and for n > 3 


B 


P ^{C P + C P - P B } 
0 n 0 n-1 1 1 n-1-^ 


CB16d) 


oo 

For sufficiently large |t| the expansion ® ^ con- 


h=0 


verges [10, pp . 54,55]. It is not important for the 
present analysis to determine how large |x| must be to 
ensure convergence. It should be noted that the trans- 
formations E, F, diagonal ^1, 1, 1, , P(x), and P(x) 

are nonsingular for X ^ X , — - k 0 , and |x| 

12 Y 

sufficiently large. 

The differential equation CB15) can be used to 
develop formal expansions about 5 = “ • The matrices 
Bq and B^ suggest attempting solutions of the form 
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L -(n+2A^) 
n=0 ^ 


(B17a) 




L - ( n+ 2X ) 

Ea T 2 

n-O 


(B17b) 




2QT 

L — — n / — 

V r t" e-^ 
n=0 ^ 


(B17c) 


and 


Z ^ (t ) 
-4 , L ^ 


2ax 


2 “'" /7 

S T e ^ 


n=0 


(B17d) 


where L = «> . For a finite value of L the formal 
truncated solutions are obtained. NOTE: The circumflex 

^ /s 

( ) is introduced to denote that ^(x) is only a formal 
expansion, that is, the full infinite expansions may 
diverge. A vector ^(x) without the circumflex will 
denote an actual solution of (BIS) . 

If formal solutions (B17) are substituted into the 
differential equation (BIS), then the following relations 
are obtained: 


1 

0 

0 

0 


(B18a) 



- 174 - 



(B18b) 


CB18c) 



If the constants in (B18) are all set equal to one, 

then the lead asymptotic vectors, , r ^ , and , 

recursively determine the remaining vectors, 

r , and s , for n = 1,2,... . 

— n ’ — n 

If 
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l^n 

2‘^n 

3^n 

4*^n 


(B19b) 


-0 


0 

0 

1 

0 
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2^n 


3^n 


4’^n 


(B19c) 
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0 

0 
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J 


and 


-n 


l^n 


2^n 


3^n 


S 

4 n 


(B19d) 


then 


iPn 


1 n- 1 
n 


(B20a) 


2 n- 1 


2 Pn 2^2 - (2A^ + n) 


CB20b) 
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3P„+1 = ^ ja^n-l * ( 2 ^ + " * l)3P„j ■ •:®20c) 

4Pn+l ■ ^ |4%-1 " (2\ " " " |)4P„| 

where 


7T 
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1 n-1 
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off -I 
3 n-1 


4”^n-l 


®2£„-1 " 


" ®„+i£o 


CBZOe) 


1 n-1 

l^n 2 X^ - (2X2 + n) ’ 

2^n-l 

2 *^n n ’ 

3*^n " 2a {3Pn-l (^^2 ^ ^ ^ ’ 

4P„ “ S j4P„-l * (^2 * n l)4'l„| ’ 


CB21a) 


(B21b) 


CB21c) 


(B21d) 


where 
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CB21e) 
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(B22a) 
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4^ I4^n-1 " ^ " 4^n ’ 


(B22d) 


where 
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2— n-1 n+1— 0 


(B22e) 
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(B23a) 
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(B23b) 
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(B23c) 
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CB23d) 


where 
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3^ 
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n- 


n- 

n- 
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1 

1 

1 

1 


B^s , + 
2— n-1 
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n + 


1-0 


CB23e) 


The formal truncated solutions satisfy the dif- 
ferential equation (B15) approximately. More precisely, 


5t 


E 

in=0 


B T 

n 


+ 0 


1 


L+1 1 




(B24) 


The recursion relations (B20) through CB23) are 
well defined if 2(X^ - 7^ integer . In addition the 

transformations E, F, diagonal |l, 1, 1, i)? PCt), and 
P(t) require that — - ^ 0 and |x| be sufficiently 

large in order for the differential equation (B15) to have 
a convergent expansion for its coefficient matrix. 
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CASE 2; The roots o£ the dispersion relation (13) 
are repeated. This is equivalent to requiring a = ^ 
in relation (13). The transformations E, diagonal 
|l , 1, 1, F, and P(x) developed for case 1 can be 
used for this case also. However, the transformation 


P(t) must be 

modified since 

II 

II 

1 

2 • 



Let 
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0 T ’ 






P(t) 

— 





(B25a) 
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(B25b) 
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, CB25c) 
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and 
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0 0 


0 0 ^ 2£_LJ. . I 

4o\ Y 2, 
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/yll 

4a \2 


2a^ + 2 ' 
Y ) 


CB25d) 


then 


where 


X(t) = PCt)Z(t) 


CB26) 


dZ(T) 

dx 
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diagonallo, 0, — , - — 
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CB28a) 
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CB28b) 
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" Pi’ ’ 


(B28c) 


and for n > 3 


n 


P^^{C + c P - P B ^ } . 

0 n 0 n-1 1 1 n-1 


(B28d) 


The rather simple form of and B^ suggests 

attempting formal solutions 


Cl. V - (n+i) 


(B29a) 


n=0 




+ (£n t)Z, ^(t) , CB29b) 


n=0 


-1 
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E 
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CB29c) 
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Z. (t) 

4 f oo ^ 


2ar 


E 

n=0 


n 


2 “ /y 
ST e ' 
— n 


(B29d) 


If ot = ^ and expansions (B29) are formally sub- 
stituted into (B27) , then the following relations are 
obtained : 

r • 

1 




(B30a) 
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^ c * 

0 

2 

1 — 
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O 

1 



9 


and 



CB30b) 


(B30c) 


CB30d) 


I£ the scalars in (B30) are set equal to one 

and ^2 ~ ^ ’ then the following recursion relations are 
obtained : 


2Pn 
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n 
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CB31b) 
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where 


I 1 

n- 1 


1 n- 1 


2 n- 1 


3 n-1 


4 n-1 


B^p _ + ••• + B , (B31e) 

2i-n-l n+1^0 ’ ^ 


2^n 


2^n-I ~ 2^n 

n 


CB32a) 
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l^n-l 2^n ~ l^n 
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(B32b) 
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(B32c) 
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CB32d) 
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CB32e) 
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CB33b) 
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CB33c) 
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CB33d) 


where 
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CB34a) 
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(B34c) 


4^n 
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n 


CB34d) 
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and 




n-1 


4 ^ 


n-1 


n-1 


n-1 


n-1 


B„s + 
2— n - 1 


+ B s „ 
n+1— 0 


(B34e) 


CASE 3: 2(X^ - X^) = integer 

If X^ and X^ differ by an odd multiple of one- 
half, then method II can be used. If X^ and X^ 
differ by an integer, then both method I and method II 
are difficult procedures to implement. 

If method I is attempted, then additional shearing 
and similarity transformations must be developed. This 
approach is considered in [8, Chapter 4]. 

The main difficulty is that the recursion rela- 
tions (B20) and (B21) are not both well defined if X^ 
and X^ differ by an integral multiple of one-half. In 
addition to the transformations E, diagonal |l, 1, 1, 

F, P(t), and P(t) , it is necessary to introduce a 
product of shearing, Sibuya, and similarity transforma- 
tions . 

Consider 


H(t) 


7l {Sim. X Sib.(r) x diagonal I-, 1, 1, 1 

i^l 1 1 \ 


(B35) 
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and ^(t) = H(t)Z_(t) . (B36) 

The shearing transformation, diagonal 1, 1, 1 j , does 

not affect in (B15) . However B^ is replaced by 

1 1 

a matrix with eigenvalues -2X^+1, "2X2, The 

matrix Sib.fT) is constructed so that the new B^ matrix 
is block diagonal with precisely two 2x2 blocks on the 
diagonal. The matrix Sim^ reduces the (1,1) block to 
Jordan form. This process is repeated N times, where 

N = 2(X^ - X^) . (B37) 


The net result is that 


dZ(x) 

dx 



(B38a) 


B, 


B, 


(B38b) 


and 

B^ has eigenvalues -2X^+N, - 2 X 2 , j, and 
— and B^ is in Jordan canonical form. CB38c) 

This is essentially the same problem as case 2 since 

- 2 X 2 = -2X^ + N . (B39) 

This procedure of developing several shearing, 


Sibuya, and similarity transformations is necessary only 
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i£ full infinite formal expansions are desired. However, 
for the purposes of this research it is necessary only to 
be assured that the lead terms in the asymptotic expan- 
sions, with algebraic growth, correspond to inviscid 
solutions and that there exist actual solutions of (25) 
which agree with the lead asymptotic terms to arbitrarily 
many significant figures as r along arg t = . 

If the transformation H(t) defined in (B35) is omitted, 
then the only loss is that full infinite formal expan- 
sions are not obtained because the recursion relations 
(B20) and (B21) are not both well defined. Since no com- 
putations for this case are necessary (see comments in 
Section 5.1), it is not worthwhile to pursue this case 
any further. 

The development of formal expansions about the 
irregular singularity 5 = 00 appears to be quite com- 
plex. A direct approach (method II) which avoids the 
complicated transformations of the differential equation 
(25) will now be developed. The success of this approach 
is based on the observation that it is necessary only to 
determine multiples of the boundary layer and transition 
layer solutions, in equations (32) and (33), in order to 
solve the viscous problem. 

Consider the differential equation (25) which can 
be written in the form 
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diCO 

d? 


ACOyCn 


CB40a) 


where 


A(C) 


II 


Y 3 


Y - 1 / n 


t, k(? - 1/3) 


/ kg \ J / g - l/3 \ 

— U - 473; - 4/3; 


g + k^ 
Y 

g - 4/3 


-(r^TT 


CB40b) 


The characteristic equation of A(|) for large ? has the 
following approximate roots : 



^1 

X 

(B41a) 



5 ■ 

(B41b) 



~ j 

(B41c) 

and 

0 

/yC 

~ 0 

~ > 

(B41d) 

where 

4 

and X^ are 

/yf 

the roots of 

the dispersion rela- 


tion (13). In addition, the eigenvectors corresponding 
to the roots Z. tend toward constants vectors as 

X 

5 ^ . The eigenvectors corresponding to and 

tend towards linearly independent vectors if the roots 
of the dispersion relation (13) are distinct. The 
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eigenvectors corresponding to and tend towards 

linearly dependent vectors. Relations (B41) suggest 

2 5 

attempting solutions o£ the form 


/\ 





L 

n=0 


- (n+X . ) 
1 


•n 


(B42) 


for i = 1 or 2 and L = °° . If formal solutions with 

/ ±2a^\ 

exponential growth le are substituted into (B40) , 

then the recursion relations are not well defined. Thus, 
method II does not yield formal solutions with exponen- 
tial growth. However, there is no difficulty in deter- 
mining formal solutions of the form (B42) . 

It will be convenient to substitute directly into 
(22) rather than (B40) . Consider 


U 


W, 


n, 1 


- (n + X , ) 
1 


n=0 


(B43a) 


where 


n , 1 


1 n , 1 

■ 

2 n , X 


(B43b) 


2 5 

The double circumflex (^)in relation (B42) denotes 

/s 

that ^ formal solution of (25) which is 

obtained by direct substitution into the original dif- 
ferential equation (method II) . 
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and 


a 

— n , 1 


A ■ 

1 n , 1 


- (n + X . ) ^ i- 
^ 1 1 n , X 

• 

2 n , X 

- (n + X . ) 

x"^ 2 n , X 


(B43c) 


Substituting CB43) into (22) leads to 



(B44a) 


or 



X . 

X 


- k" - 1 id) 

T Y 


(B44b) 


and 



1 

0 , X 

2 / 1 2 
o /y - K 


- ?) 


CB4 5a) 
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o 


Y 


k ( X ^ + n + 1 


Y - 1 
Y 


kX. +n+ 1 - — 
1 Y 


(X^ + n + 1 ) (X^ + n) + Y 




n+1 , i 


(X. + n)2 





4(X^ + n)^ 


I . 

n , 1 


(B45b) 


The matrix coefficient of £ , , . is nonsingular 

in (B45b) for all n if and only if X^ and X^ do not 
differ by an integer. Thus, it is possible to determine 
formal solutions of [22) via the recursion relation (B45) 
for ■ ^2 ^ integer . 

CASE 4: ^ - k2 = 0 

Y 


Different normalizations of the lead vectors, I . 

u , X 

[i = 1 or 2) , are required. 




(B46a) 


and 




X^ = ^ , (B46b) 


and the remaining vectors (5,^ ^ , i = 1 or 2) can be 
determined from the recursion relation (B45b). 
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Tlius , there is little difficulty in determining 
formal solutions of the differential equation (22) about 
the irregular singularity C = - However a fundamental 

set of formal solutions has not been obtained. Approxi- 
mate multiples of the boundary layer and transition layer 
solutions can be found instead of determining two addi- 
tional formal expansions. 

Suppose, due to (B41) , it is assumed that there 
exist two solutions, y 3 (?) and (C) , of the differential 
equation (25) which satisfy 



MisCali 



(B47a) 


and 



-c 


3 

e 


-Re2a 



< 


1/4(5) I 



(B47b) 


for some positive constants M^, M^, M^, , c^, c^, c^, 

and c' . 

4 

If such solutions exist, then there should be 
very little difficulty in determining multiples of these 
solutions numerically. The solution y_ 3 (?) grows more 
rapidly than }^(5) for increasing |^| and arg ^ = -j . 
In addition, there are solutions y^ 3 ^(?) and which 
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A A 

are asymptotic to I-2 00^^^ hence grow in 

-ReX^ ' -ReX 

norm essentially like |Cl ^ and |^| respectively. 

Thus, y^^^^ grows more rapidly than and 

1^1 increases along the ray arg 5 = • 

Consider a numerical integration in the direction 
of increasing |^|. For an arbitrary initial vector it is 
expected that the initial vector will contain nonzero 
multiples of all four solutions, y^CC), ^nd 

Since 73(0 dominates in growth as |c! is 
increased, a numerical integration of ( 25 ) over a suffi- 
ciently large |^| interval in the direction of increasing 
1^1 should result in nearly a multiple of ’ Simi- 

larly, if a numerical integration of ( 25 ) is performed in 
the direc ion of decreasing |C|, then a multiple of y^(C) 
can be obtained since y^(?) has dominant .growth for 
decreasing |^| and arg ^ . 

The procedure (method II) leads to very good 
agreement with the asymptotic solutions obtained by 
method I, that is, multiples of 73(0 and y)^(?) are 
obtained. The agreement of these two approaches provides 
a good check on all the transformations performed in 


method I . 
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B.2 PROPERTIES OF THE ASYMPTOTIC SOLUTIONS 


In this section the formal asymptotic solutions 
developed in Section B.l are investigated. All of the 
various properties of the asymptotic expansions mentioned 
in Sections 3, 4, and 5 are developed. 

Associated with the formal asymptotic solutions in 


Section B.l there are four distinct rates of growth, 


26 


20T 


20T 


•2X- 


•2X 


2. ^ 


— ^ /\y ^ vy j z. 

, T ^ y T e , and t e 


/y 


J 2a-J^ i - 20 ^ 

•4^ lY Ir 


-^1 "^2 

5 \ re and re 


or equivalently 


Recall that 


lea 


and T = /^ ; thus 


-2X. 


Xi X^z 

(iea) e 


(B48a) 


-2X. 


(lea) e 


(B48b) 


1 20 

2 /y 

T e = 


■ z \ 4 


iieo 


,2a Iq 
exp ^ — 


•z/2^ 


/y Viee 


(B48c) 


2^ -20 

j 2 /y 

and T e ' = 


■z \4 


U CO 


exp 


2o /e 


■z/2\ 


/y Vico 


(B48d) 


It is tacitly assumed that X^^ X^ and 

o /y - k / 0 . The several special cases will also be 
investigated . 
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1 

In Section 3.2 the scalar growths x^e were 
investigated and it was shown that the solutions with 
these asymptotic growths are important only in the 
boundary layer and transition layer. The scalar growths 

(ieo) ■^e and (ieo) are obviously the same as the 

scalar growths associated with the inviscid solutions. 

\ ^2 

The factors (ieo) and (ieo) are merely scaling con- 
stants for a prescribed e > 0 . In addition to the 
scalar growths it is also important to investigate the 
lead vectors in the formal asymptotic expansions. 

Let 

Tq = E diagonal(l, 1, 1, 0 )FP^Pq CB49a) 

T^ = E{diagonal (1 , 1, 1, 0) [FP^P^ + 

+ diagonal(0, 0, 0, 1 )FP^Pq} (B49b) 

and 

T(x) = E diagonal ^1, 1, 1, ^jFP(T)P(x) , CB49c) 

where the matrices E, diagonal (l, 1, 1, — j, F, PC't)? 
PCt), Pq, P^, Pq, P^ are defined in Section B.l. 
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Consider 


(O = T(t)Z Ct) (i = l,2,3, or 4) (B50a) 

X / J-J X ; Xj 

and 


= T(t)ZJt) 


CBBOb) 


where ^ = t and Z. ^ is defined in relation (^17). 
Thus, £. ^C?) is a formal solution of (25) since Z^. (t) 

is a formal solution of (B15) . If ^^(x) is an actual 
solution of the differential equation (B15) , then 
is an actual solution of the differential equation (25) . 
The formal truncated solution (?) is of the 

form 




L, 

X 

n=0 


b 

■n , 1 


-n/2 




■X . 

X 


and 




for i=l or 2 


n=0 


L+l\\ / 1 d . 20 


5 e 


3 ly 


4 


(B51a) 


(B51b) 


for j = 3 or 4 and d^ = +1 and d^ = -1 . The vec- 
tors b . and c . are determined by the transforma- 

-n , X -n , 3 

tion T(t) and the formal truncated expansions 



! tsi> 
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. ^ (t) (B17) . In particular, the vectors b . and 

i ^ L — 0 / i 

. are given by 

- 0,3 
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and 


- 0,4 






a 


0 

1 

0 

0 


(B52d) 


Notice that the lead vectors c^ ^ and c^ ^ are mul- 
tiples of each other. This partially explains the 
difficulty in determining the formal asymptotic solutions 
with exponential growth. 

The formal solutions (B42) can be determined from 

relations (B43) and (B45) . The lead asymptotic vector 

a^ . is given by 

— 0 , 1 

1 

-X , 

X 


-0 , i 



CBS 3) 


for i = 1 or 2 . Thus, the formal solutions (^42) have 

"^1 -^2 

the lead asymptotic terms Uq and 2 ^ .It has 

-Xi 

already been established that the scalar growths 5 
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-^2 

and 5 correspond to the inviscid growths e 

and e , respectively. Comparison of relations (B53) 
and (15) establishes that the lead asymptotic vectors 
a and a^ correspond to the inviscid vectors. 

Note that the first and third components of a , ? ^ 

U r X 

should be compared with U^(z) and W^(z) in (15), respec- 
tively, due to relations (24) . The second and fourth 

-X. dU.(z) 

components of a should be compared with - — t— 

dW, (z) 


and -- 


dz 


These comparisons are easily made since the 


normalizations of the lead asymptotic vectors a and 

U f 1 

a consisted of specifying the first components of 

\j f ^ 

each of these vectors to be one. The formal solutions 
(B17) are normalized by specifying the constants c^ in 
relations (B18) to be identically one for i = 1, 2,3, 
and 4 . The lead asymptotic vectors ^ ^ given in 
(B52) have a rather complicated form. However, it can be 
verified that 


Y 


-0,i 


Y 


tr O - 


CB54) 


for i = 1 or 2 . 

Thus, the formal solutions (B17) for i = 1 or 2 
are asymptotic developments of actual solutions of (25) 
and the lead terms in these formal solutions correspond 
to multiples of the inviscid solutions (15) . 
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2 2 

For the case a = yk , the normalization (B46) 
leads to a development o£ which corresponds to the 

inviscid solutions (l6) . 

One additional case must be considered, namely the 
case ~ ^2 ~ Y formal truncated solutions (B51) 

are of the form 



CB55a) 



where j = 3 or 4 , d^ = +1 , and d^ = -1 , 

Once again it is of interest to investigate the 
lead asymptotic terms which do not exhibit exponential 


growth or decay. Consider 





1 

2 



and 
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1 

1 

2 



1 

(iea) ^ 




If the first and third components are compared with 


the inviscid solution 


Ul(z) 
j^Wi (z)j 


in equation (18) , then 


~ 2 

it is easily seen that b p, , C is a multiple of this 

U , X 

inviscid solution. Now consider 



) -j 

(ta c)b„_^|c 


(iea) 


z/2|^ 


-0 , 2 


(iea)^e^^^^i £n 




(B56) 


(B57) 


where 
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- 0 , 2 


Y ( 4 - Y 
Y - 1 \ 4 

T 


4Cy - 1) 




Y 


4 (y - 1) 


'XUk2 - 2- 
Y 


(B58) 


The first and third components of ^ „ do not compare 


with the expression for 


UsCz) 


in relation (18) . How- 


ever, the two-dimensional vector which is constructed 

from the first and third components of 
1 

z 




2 bp is a solution of the differen- 


tial equation (12) . Hence 

\ 1 

1 I ~ 2 

1—0 2^2' ^^—0 l(^ compares with a multiple of an 

inviscid solution. 

Thus, in all cases the lead terms in the formal 
expansions with algebraic growth correspond to multiples 
of the solutions to the inviscid differential equation 
(12), that is, the lead terms are merely inviscid solu- 
tions . 
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In order to distinguish the asymptotic expansions 
for the boundary and transition layer solutions, it is 
necessary to determine additional terms in the asymptotic 
expansions. The relations (B22) , (B23) , (B49) , and (B51) 
yield 




- 1,3 



and 




- 1,4 


- 1,3 


CB59b) 


where p^^ is the (4,3) element of defined in rela- 

tion (B14) . The same result (B59) is obtained for the 

case ^1 ~ ^2 ~ 2" ■ 

LEMMA B1 : The vectors b^ b^ c^ and c^ ^ 

defined in relations (51) are linearly independent if 

2 2 

~ ^2 ^ integer and a /y - k / 0 . 
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PROOF 


det £^^3. £3^3] 


_ /? 


Y 


a \y - 1 


)%X3 - X3>|r 


^ 2 \ 2 
g Y 

Y ^ 


a 


CB60) 


If -0 , 1 ’ - 0 , 2 ’ -0 , 3 ’ -1 , 3 ^ inear ly dependent, 

then CB60) implies 


.2 0 ‘ 


(B61a) 


or 


CB61b) 


or a = 0 => = 1 , = 0 . (B61c) 

The hypothesis of lemma B1 excludes relations (B61) . 

Hence, b^ c^ and c^ ^ are linearly inde- 

pendent. Q.E.D. 

For computational purposes it is important to 
avoid near-linear dependence or ill -conditioning of the 
vectors b^ b^ c^ and c^ ^ . Although lemma B1 
implies that the vectors b b c , and c are 

linearly independent for all k , it was found that for 
k > 5 , near-linear dependence destroyed the calcula- 
tions. The error check, step 5 in the algorithm developed 
in section 4.1, implied that no significant figures were 
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obtained for these calculations (k > 5) . To avoid this 
difficulty the range of k values was restricted 
(k < 1.5). Similarly, for ^ ^2 ^ \ Lamb 

o p 

wave (a /y ■ k =0), good accuracy was obtained for 
k < 1.5 , that is, the problem of ill -conditioning was 
avoided . 

Relations (B52) and (B59) imply 


and 


BLSOL(g) « 


- 0 , 3 


+ 




TLSOL(g) 



+ 




(B62) 


(B63) 


for large ? and arg ^ . Notice that the form of 

the vectors c and c implies that the modulus of 

the third component of the boundary layer solution is 
much smaller than the modulus of the first component. In 

TT 

fact, as ^ 00 along the ray arg ^ ■ 


third component of BLSOL (C) 
first component of BLS0L (XJ 

For e = 10~^^ the ratio (B64) was approximately 10 ^ at 
z = 0 . Thus, in most cases the boundary layer solution 
which arises in equation (33) affects only the 


= . (B64) 
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9 7 

horizontal component of the velocity. The boundary 
layer solution reduces U(z) to zero in the boundary layer 
and has a small effect on W(z) . If the boundary layer is 
ignored, then the constant d^ in equation (33) can be 
accurately determined so that the kinematic boundary con- 
dition W(0) =1 is satisfied, that is, it is only 
necessary to satisfy 

Zvptq) = • (b65) 

Equation (B65) will determine d^ to about five signifi- 
cant figures for e = lO"^^ . Of course only the kine- 
matic boundary condition W(0) = 1 will be satisfied by 
(B65) . 

Thus, the boundary layer solution has negligible 
effect on W(z) for all z , including z = 0 . In 
addition, it was shown in Section 3.2 that the boundary 
layer solution has negligible effect on U(z) outside a 
thin boundary layer near z = 0 . 

In Section B.l two distinct methods were developed 
for computing formal solutions which have lead terms cor- 
responding to multiples of the inviscid solutions. For 
computational purposes it is necessary to consider addi- 
tional terms in the formal expansions. Quite naturally 

2 7 

The only exception is the Lamb wave. 



-207- 


the question arises as to whether or not the remaining 
terms in the two distinct developments (method I and 
method II) are related in any way. 

So far it has been shown that 


c. X a 

X —0 , X 


(B66) 


where the vectors b , and a . are defined in (B51) 
and (B42) , respectively. 

THEOREM B1 : If relation (B66) is satisfied and 
2(Xi - X^) 7^ integer and a^/y - 0 , then 



for i = i or 2 and n = 0,1,2, ••• . 

NOTE: If theorem B1 can be established, then 

there is very little difference in using defined 

in (B42) or y, (5) defined in (B51) for computational 

““X f Zxj 

purposes . 

PROOF: Let be a solution of the differen- 

tial equation (25) such that y^(?) is asymptotic to 
y. ^(5) as 5 oo along the ray arg ^ = -y , that is, 

MZiCC) - Z.i,L(5)M = ^ jlcl ^ , (B68) 



-208- 


£or i = 1 or 2 and all L = 0,1,2, . For a proof of 
the existence of an actual solution of the differen- 
tial equation (25) which is asymptotic to see 

[8, Chapter 5] . Similarly, let Y^C?) be a solution of 

the differential equation (25) such that Y^(5) is asymp- 
/\ « 

totic to y. ^(?), that is, 

IIYJO = o(r Ul ^ , (B69) 

for i = 1 or 2 and all k = 0,1,2,3,-** . 

Consider w^(C) defined by 

w.CO = Y,(?) - . CB70) 

Since (25) is a linear homogenous differential equation, 
w. Cn is a solution of (25) . The hypothesis and rela- 
tions (B68), (B69) , and (B70) imply 

/ -i-\ -ReX 

I |w. (O I I = 0^5 Vl?| 

If 



0 


(B71) 


then 


W]_(?) ~ 




c X a 1 
1 - 1 , 1 ^ 




(B72) 
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I£ 2(X^ - X^) / integer, then (B71) and (B72) imply 

14(^3 are linearly 

independent a-nd y^CO ate defined in relation 

(B50) j , since each of these solutions of the differential 
equation (25) has a different asymptotic growth. How- 
ever, this is impossible since a system of four first- 
order, linear differential equations has only four 
linearly independent solutions. Relation (B72) is a con- 
sequence of relation (B71) ; thus (B71) must be invalid or 

All the coefficients in the asymptotic development 
(B72) must be zero; otherwise w^^C?), ’ 

and y_^CC) are linearly independent. Similarly, 






c X a ) 
2 -1,2'' 




(B73) 


If any of the coefficients in the asymptotic develop- 
ment (B73) are nonzero, then > 

y (^), and y. (C) are linearly independent. Since this is 
impossible it follows that 



(b 


2 , 1 


^ 2-1, d 


0 


and 



(b 


2 , 2 


S " ^1,2^ 


0 


Q.E.D. 
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Thus , for computational purposes there is little 
difference in using either the formal expansions (B17) 
and transformations (B49) to obtain approximate values of 
the solutions of (25) with inviscid growth or the expan- 
sions (B42) . 

Since the asymptotic solutions are used for com- 
putational purposes, it is necessary to obtain some error 
estimates. Consider 


Mzpu - ii,L(UII - 0(T-^-p||Z._^(T)|| , (B74) 

as T ^ , where Z, (x) is defined in (B17) and (x) 

is a solution of (B15) . If e > 0 is prescribed, then 
there corresponds a constant M(e) such that 


-L + 1 


- e 


-2X. 


L+1 




-L+1 


+ e 


L+1 


-2A. 


(B75a) 


l^L+l 


- e 


-2X, 


L + 1 




3-Ltll I * ^ 

I _ I L+1 


•2X. 


, (B75b) 
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-L+1 


L + 1 


T e 


20T 
1 Re 

/7 


and 


1 -Re 


-L+1 


L + 1 


I 2 

T e 


20T 




II 1 Re 

JIl+JI " " I ,5 

■ T e 


2aT 

/7 


L+1 


(B75c) 


Z(t) - Z (t 
-4 ^ — 4 , L ^ 


11^ 1 -Re 

S + £ — 

-L + 1 I I I I 2 

• T e 


20T 

/? 


L+1 


' ' (B75d) 

for |t1 > M(e) and arg i • It is difficult to 

determine M(e) . However, an indication, not a bound, of 
the error can be found by setting e = 0 in (B75). For 
a specific value of t it would be reasonable to trun- 
cate the expansions (B17) at a value of k such that 
r(x,k), defined by 


r(x,k) 


I IRv 


max 


llaJI 


id I IliJI 


(B76) 


is a minimum; that is, the expansions (B17) should be 
truncated at L = N , where 


foi" all k = 0,1,2,*»* (B77) 

and X is held fixed. Of course N will vary with the 
value of X prescribed. 
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Notice that r(x,N) provides an estimate o£ the 
error. For a specific value of x there will correspond 
a certain degree of accuracy. However, it is reasonable 
to require a certain accuracy from the expansions (B17) 
rather than accept any accuracy that is possible at a 
prescribed value of x . For example, if it is desired 
to satisfy 

r(x,N3 < lO""^ , (B78)^^ 


then X and N must be determined. Due to many calcula- 
tions it was empirically established that (B78) is satis- 
fied if 




(B79a) 


and 


N « 20 


(B79b) 


It is possible to verify 
recursion relation (B45) or 


(B79) by examining the 


£ 


n + 1 , i 


^ + 0(n) 0(n) 


I 


0(n) 


0(n) 


n , 1 


(B80) 


2 8 n • 

Most of the calculations summarized in Section 5 
were performed on an IBM 7094 in single precision (about 
eight significant figures) and relation (B78) was 
satisfied. 
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A crude estimate can be obtained for N and |x| by 

n I 2 

examining ^ Recall that N = 20 for 


expansions (B17) 
formal expansion 

(1-4) 


is equivalent to summing 10 terms of the 
(B42) . It is easily verified that 


( 10 ! ) ^ 

20 


-10 


< 10 ’ 


(B81) 


if 1^1 = (^) • CB82) 

More generally, if it is required that r(T,N) satisfy 

r(x,N) < 10"P , (B83) 

then 

A 

j X j (A is some positive constant) (B84a) 

and N varies only slightly, (B84b) 

Since the function r(x,L) does not vary signifi- 
cantly for L near the value N defined in (B77) , it is 
not critical to truncate the expansions (B17) at precisely 
the minimum term. Hence, it is possible to fix the value 
of N with very little change in the determination of 
the initial vectors for the canonical numerical integra- 
tion procedure. 
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In Section 3.2, equation (32) was developed. 
Suppose the same notation is used for the solutions of 
the differential equation (25), namely DC^ (g) , DC2(5), 
INV ^ (g) , INV ^ (g) , and TLSOL (g) , and suppose the scalars 
®2’ defined in (32). It is of 

interest to investigate the solutions of (32) as e 0 . 
Thus, it is natural to investigate the dependence of e^, 
®2’ ®3’ parameter e for fixed values of 

a and k . 

THEOREM B2: If the solutions of the differential 

equation (25), which are asymptotic to multiples of the 
inviscid solutions, are required to satisfy 


■n’ , , . 



INV^(g) 

~ e^ T(t)Z^^^(t) 

(B85a) 

and 

INV2(C) 

7T / -V V 

yIm(A ) 

~ e T(t)Z2^^(t) , 

(B85b) 

where 

and A ^ 

are the roots of the dispersion 

rela- 

tion (13) , 

T(t) is 

defined in (B49) , and Z, ^(A) is 


defined in 

(B17), 

then the scalars e, , e , and e 

1 ’ 2 ’ 3 

are 

invariant 

as e ^ 

0 if DC^(g), DC 2 (?), INV^(g), and 


TLSOL (g) are linearly Independent. 

PROOF: The vectors DC^ (g) and DC^ (g) are uniquely 

defined. The vectors which satisfy (B85a) form a one- 
parameter family of solutions such that any two members 
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o£ the family differ by a multiple of the transition 
layer solution. Similarly, the vectors which satisfy 
(B 85 b) form a one-parameter family of solutions. Thus, 
consider 

e^^iCC) + 62^2^0 + + e^ TLSOL (g) = 1^2(5) 

and 

f^DC^CO + ^ f 3 {IW^( 5 ) + A TLSOL tg)} + f JLSOL (g) 

= INV2 (f ) + B TLSOL (g) . (B86) 

The theorem is established if for arbitrary scalars A 
and B , ^ ~ 1 > 2 , and 3 . However, equa- 

tion (B86) can be rewritten as 

- £2^2^^^ " £ 3 lNV ^(0 

+ {f^ + f3A - B} TLSOL (g) = INV2(g) 

If DC^C?), DC2(?), INV^(g) , and TLSOL (g) are linearly 
independent, then 13^ = 6^, £2 ~ ^3 ~ 

+ fjA - B = . Q.E.D. 

In general, if conditions (B 85 ) are satisfied, 
then the scalar e^ can vary. In addition, if INV ^ Cg) 
and INV 2 (g) are reduced to canonical form at ^ ~ J > 
then the scalar e^ is also uniquely determined. 
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THEOREM B3: 
I£ 


DC^C?) 01" DC^ eg) is asymptotic to a nonzero 
multiple of the boundary layer solution 


(B87a) 


and 


whenever C?) +<^ 2^2 asymptotic 

to a nonzero multiple of the boundary layer 
solution, then 

where and f^ are nonzero if 

Ic^l + IC 2 I > 0 , 


(B87b) 


then the viscous problem has one and only one solution 
corresponding to every sufficiently small e > 0 for 
the following cases 


a . 




7^ 0 


b . 



1 

2 


29 , 

C. A2 < 



Y = 1 .4 , 


^^For case c, hypotheses (B87a) and (B87b) are modi- 
fied slightly. 
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where and X^ are the roots o£ the dispersion rela- 

tion (13) and y_. ^(5) is defined in (B50) . 

PROOF; Due to hypothesis (B87a) , either DC^ (g) or 
is asymptotic to a multiple of the boundary layer 
solution and, hence, simply by interchanging subscripts 
it is always possible to assume DC^ (g) is asymptotic to a 
multiple of the boundary layer solution. Thus, without 
loss in generality assume DC^ (g) is asymptotic to a mul- 
tiple of the boundary layer solution, that is, 

where a ^ 0 and y defined in (B50) . 

Consider 

For a unique value of b it is possible to eliminate the 
boundary layer solution from DC^ (g) . If DC^ (g) is not 
asymptotic to a nonzero multiple of the boundary layer 
solution, then b = 0 and if DC^ (g) is asymptotic to, say, 
^^Z -3 00^5^’ then b = -^ . Thus, for the proper choice of 
b , asymptotic to a nonzero multiple of 

the boundary layer solution and hypothesis (B87b) implies 

where c and d are nonzero. 
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Now let and dc^CC) be constructed from the 

first and third components of DC (5) and DC (Q , respec- 

X ^ 

tively. Clearly (see theorem 1 in Section 3.2), if it 

• — «"S./ 

can be shown that Hc’^(5) and dc^CO are linearly indepen- 
dent for sufficiently large ^ , then the theorem is 
established. 

2 

CASE 1: ^ q 

1 2 2 ’ Y 

Relations (B52) , (B59) , (B88) , and (B89) imply 



Hence , 

det [dc^(C), (Tc^COl 




“•4 

e ' 



+ 0 ( 1)5 


(B91) 
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1 a 2 

and by assumption ^2^2’ ~ ~ ^ 0 > e ^ 0 3-n-d- g ^ 0 . 

Thus, relation (B91) implies that for sufficiently large 
5 , dc^C5) and dc^C?) are linearly independent. 

CASE 2: = ^2 " 7 

Relations (B55) , (B56) , and (B89) imply 

(£n 5)5 ^ . 


dc^ (5) 


2 / 

a /y 


o (1) 


kd - i 

2 Y 


0 ( 1 ) 


det [dc^(5), cTc^CC)] = 


Hence, relation (B92) implies that 3^^ (5) and dc^CC) ate 
linearly independent for sufficiently large 5 • 

CASE 3: 0^/y - = 0 , y = lA 

For this case the formal expansions were obtained 
by method 11. However, by investigating the eigenvectors 
of the matrix A(5) defined in relation (B40) , it is 
expected that relation (B90) is still satisfied, that is, 
M^(5) satisfies (B90) . Thus, hypothesis (B87a) is 
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slightly modified and relation (B90) is assumed to hold. 
In addition, hypothesis (B87b) must be modified; that is, 
assume 

~ ’ CB93) 

if c^DC^ (g) +c^DC^ (g) is not asymptotic to a nonzero mul- 
tiple of the boundary layer solution, where y, is 

defined in relations (B42) , (B43) , (B45) , and (B46) . 

Thus , 




and 

det [dc^(0, ^c^C?)] 


However, e ^ 0 (relation (B90) assumed), m ^ 0 , 
and ^ ^2 ^ ^ T (recall y = 1>4); thus, 

dc^(C) and are linearly independent for suffi- 

ciently large 5 . Q.E.D. 
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THEOREM B4: 
If 


and are complex roots of the 

dispersion relation (13) , 


(B94a) 


DC (g) or M (5) is asymptotic to a nonzero 

^ (B94b) 

multiple of the boundary layer solution, 


Y_ equals a multiple of the inviscid 

1 , (J 

solution with upward energy flux and y 

(B94c) 

equals a multiple of the inviscid solution 
with downward energy flux, 

and 


whenever c^DC^ (g) +c^DC^ (g) is not asymptotic 
to a nonzero multiple of the boundary layer 
solution, then 


c^^^(0 + c^DC^CO ~ a y^^^COe 




+ K y^ (^)e 
R^2 


jImU2) 


(B94d) 


where a ^ 0 if |c I + |c | > 0 , jx | < 1 and 

LA R 

Y_. is defined in (B50) , then for sufficiently small 

e > 0 the viscous problem has one and only one solution. 
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PROOF : Without loss in generality, assume DC^ (g) 

is asymptotic to a nonzero multiple o£ the boundary layer 
solution; that is, 

~ 43,oo(?) ^ CB95) 

where y; (g) is defined in CB50) and c ^ 0 . 

For a unique value of b , the vector DC ^ ( g ) 
defined by 

satisfies 

^Im (X ) ^ yim (X ) 

where a 0 and [k | < 1 (due to hypothesis (B94d)). 
Let dc^(g) and dc^CC) consist of the first and third com- 
ponents of and > respectively. Relations 

(B52), (B59) , (B95), and (B96) imply 



dc^(g) 




? 
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and 


<ri,C5) 



1 


+ 


K 


R 


Y 

¥ 


(1 


0 


^2^ 



■^2 

^ e 



where c f 0, a ^ 0 and |k | < 1 . 

R 


NOTE : 


1 - = 11-^2! ^ ° 


and 


? e 


-X |lm(X2) 
^ e 


Ul 


1 

2 


Im (X^) 
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Thus , 


det [dc^ (5) , dc^ (5) ] 


Y 


(y - l)a^k 


Y 


X 111 




X 



0(1) U I 



(B99) 


For |k I <1 , relations (B97) , (B98) , and (B99) imply 
that and dc^CC) are linearly independent for suffi- 

ciently large C • Hence, the viscous problem has a 
unique solution for sufficiently small e > 0 . Q.E.D. 

It has been tacitly assumed that k > 0 . However, 
the limiting case k = 0 is of value. If a is bounded 
away from zero as k -> 0 , then relation (B53) implies 
that the inviscid solutions tend towards solutions with 
only vertical motion. Thus, as k ■> 0 the wave motion 
approaches the more restrictive case of vertical oscil- 
lations of an isothermal atmosphere. The limiting case 
k = 0 has already been solved by Yanowitch [3]. 
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